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ABSTRACT 

The  analysis  of  vibration  generation  and  transmission  in  shipboard 
machinery  installations  ha£  defied  rigorous  mathematical  description. 
After  certain  simplifying  assumptions  have  been  made,  the  vibration  response 
of  the  installation  can  be  predicted  analytically.  Theeobject  of  this  thesis 
is  to  demonstrate  the  method  by  which  the  simplified  system  can  be  analyzed 
to  give  vibration  trends  and  characteristics.  The  method  incorporates  the 
four  pole  parameter  technique  of  analysis  and  a  systematic  study  of  the 
effects  of  installation  modifications. 

The  four  pole  parameters  of  a  system  are  those  expressions  which 
inter-relate  two  basic  input  and  output  quantities.  For  vibrating  mechanical 
systems  these  quantities  are  force  and  velocity.  The  four  pole  parameters 
must  be  derived  or  measured  for  the  mechanical  elements  which  transmit  the 
force  and  velocity  through  the  installation.   In  a  shipboard  installation 
these  elements  are  the  machine,  mounting,  foundation,  hull,  and  water. 
After  the  four  pole  parameters  are  obtained  for  these  five  elements, 
matrix  notation  and  combination  permit  a  transmission  study  of  the  entire 
five  element  system.   The  computational  assistance  of  a  digital  computer 
is  an  indispensable  aid  in  this  portion  of  the  study. 

Numerous  conclusions  and  recommendations  were  obtained  from  the 
analysis  and  from  the  results  of  the  systematic  study.  The  most  important 
of  these  were: 

(1)  For  the  purpose  of  vibration  analysis  a  complex  machinery 
installation  can  be  represented  by  a  relatively  simplified 
system.  This  representation  and  the  use  of  the  four  pole 
parameter  technique  are  sufficient  to  establish  trends  and 
characteristics . 

(2)  The  stiffness  of  the  resilient  mounting  has  a  dominant  effect 
on  the  vibration  transmission  at  all  frequencies. 

ii 


(3)  A  Portsmouth,  mounting  with  decreased  stiffness  and  increased 
damping  and  synthetic  rubber  damping  applied  to  the  foundation 
were  the  most  effective  installation  modifications  in  providing 
overall  transmission  reduction. 

(k)     As  new  ideas  for  reducing  vibration  transmission  are  intro- 
duced, they  should  be  investigated  by  means  of  the  four  pole 
parameter  technique. 

(5)  A  model  of  the  approximate  system  of  this  thesis  should  be 
constructed  to  permit  experimental  verification  of  the  four 
pole  parameter  technique. 


Thesis  Supervisor:  John  R.  Baylis 

Title:  Associate  Professor  of  Naval 
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CHAPTER  I 


INTRODUCTION 


1.1  Importance  of  Vibration  Transmission 

The  problem  of  noise  control  is  one  of  major  importance  to  naval 
engineers.  A  ship  with  a  high  noise  level  operates  under  the  joint 
handicaps  of  being  easily  detectable  by  the  enemy's  listening  devices  and 
of  reducing  the  effectiveness  of  its  own  sonar.  Noise  control  implies 
the  control  of  vibrations  since  acoustic  or  noise  energy  in  a  medium  is 
synonymous  with  vibration  transmission  in  that  medium.  A  thorough  under- 
standing of  the  vibration  transmission  characteristics  of  shipboard 
installations  can  provide  the  solution  of  the  noise  control  problem. 

1.2  Current  Developments 

Spurred  by  the  Navy's  need  for  quieter  submarines,  many  facilities 
and  laboratories  have  been  devoting  time  and  effort  towards  finding  a 
solution  of  the  noise  problem.  Many  studies  are  currently  in  progress 
seeking  mechanical  means  of  reducing  vibrations  In  machinery  installations . 
Much  of  this  is  being  financed  by  Navy  contracts;   thus,  most  of  the  major 
new  developments  have  shipboard  applications . 

The  development  receiving  the  most  application  at  present  in  the 
Navy  is  the  rubber-type  resilient  mount.  These  Navy  mountings  have  been 
developed  by  the  Engineering  Experiment  Station,  Annapolis,  Maryland  and 
the  Portsmouth  Naval  Shipyard,  Portsmouth,  New  Hampshire  (15)-  The 
effectiveness  of  resilient  mountings  has  been  studied  by  Mr.  A.  0.  Sykes 
of  the  David  Taylor  Model  Basin  (13,  lk) .  Studies  have  also  been  made  of 
the  properties  of  other  rubber-like  materials  and  of  their  effectiveness 
when  used  as  mountings  (17)« 

It  is  becoming  more  apparent  to  interested  engineers  that  the 
other  elements  in  the  vibration  transmission  path  from  source  to  water  must 

(1) 


(2) 
be  investigated.  Work  has  been  done  in  determining  vibration  character- 
istics in  machinery,  the  most  common  vibration  source  (12).  Machinery 
foundation  vibrations  and  techniques  of  damping  them  are  being  experimen- 
tally evaluated  by  the  Sound  and  Vibration  Section  at  the  Electric  Boat 
Company  (20) .   This  work  has  emphasized  the  importance  of  the  foundation 
in  the  vibration  transmission  path. 

The  final  item  in  the  vibration  transmission  path  is  the  ship's 
hull.  Measurements  of  the  vibration  response  of  submarine  hulls  have 
been  made.   This  empirical  data  gives  a  possible  means  of  predicting  the 
interacting  effects  of  foundation,  hull,  and  water  (ll). 

It  has  been  recognized  that  the  time  to  apply  noise  control 
measures  aboard  ship  is  in  the  preliminary  design  and  construction  stages. 
This  concept  has  been  applied  to  the  SSN  593>  Thresher,  to  be  commissioned 
this  summer.  A  consistent  set  of  specifications,  setting  allowable  noise 
and  vibration  levels,  have  been  developed  for  this  ship. 
1.3  The  System  Under  Study 

There  are  several  means  by  which  noise  or  vibration  aboard  ship 
can  be  transmitted  from  its  source  to  the  water.   If  an  operating  machine 
is  considered  as  a  typical  noise  source,  its  vibration  can  reach  the  hull 
and  water  by  several  paths.  These  are  airborne  transmission,  supporting 
structure  transmission,  and  transmission  through  connected  piping.   The 
vibration  paths  are  best  illustrated  by  the  component  block  diagram  in 
Figure  I . 


(3) 
Figure  I 
System  Block  Diagram 
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The  initial  step  in  formulating  an  analytic  approach  to  this 
transmission  problem  is  to  assume  that  each  of  these  paths  can  be  treated 
separately.   This  thesis  is  limited  to  an  analytic  study  of  the  trans- 
mission through  the  supporting  structure.  The  path  under  study,  as 
indicated  in  Figure  I,  consists  of  machine,  mounting,  foundation,  hull, 
and  water. 
1.4  The  Analytic  Approach 

The  exact  describing  equations  of  the  transmission  phenomena 
through  the  above  mentioned  mechanical  elements  would  be  an  ideal  analytic 
tool  for  the  investigation  of  this  problem.  The  actual  system,  however, 
defies  any  clear-cut  description,  and  to  obtain  the  describing  equations 
is  a  Herculean  task.  The  system  is  non-linear  and  possesses  many  degrees 
of  freedom.   The  vibration  of  most  machines  is  of  a  random  nature.   In 
addition,  these  machines  are  mounted  on  complex  structures . 

To  circumvent  the  above  complexities,  simplifying  assumptions 
must  be  made  to  obtain  a  workable  mathematical  model.   These  assumptions 
are  completely  enumerated  in  Chapter  II  of  this  thesis.  The  simplified 
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model  is  not  intended  to  give  the  complete  picture  of  vibration  trans- 
mission in  the  actual  installation. 

The  object  of  this  thesis  is  to  make  a  systematic  study  of  the 
analytically-described,  simplified  system.  It  is  hoped  that  from  this 
simplified  study  vibration  trends  and  preferred  design  practice  for  the 
complex  machinery  system  will  become  available. 


(5) 

CHAPTER  II 
PROCEDURES 

2.1  The  Preliminary  Experiment 

The  authors  first  became  interested  in  an  analytic  approach  to 
vibration  transmission  in  shipboard  installations  while  on  temporary  duty 
at  the  Admiralty  Engineering  Laboratory,  West  Drayton,  England.  A  series 
of  experiments  were  being  conducted  there  to  determine  isolation  charac- 
teristics of  Admiralty  resilient  mountings.  Results  had  been  obtained 
from  one  of  these  basic  experiments,  and  the  authors  investigated  various 
approaches  to  the  mathematical  verification  of  this  data. 

The  experimental  results  studied  were  obtained  by  driving  a  beam 
supported  by  an  Admiralty  mounting  with  a  typical  vibration  generator. 
This  experiment  is  described  in  Appendix  A-l. 

Most  previous  experiments  investigating  the  characteristics  of 
mountings  treated  the  supported  member  as  a  lumped  mass.   This  treatment 
facilitated  use  of  a  theory  involving  electrical  analogies  (14).   In  this 
case  the  beam  was  used  as  the  supported  member  because  of  its  non -rigidity. 

A  recent  paper  by  R.  D.  Cavanaugh  and  J.  E.  Ruzicka  provided  a 
method  by  which  a  vibrating  beam-spring  system  could  be  analyzed  (h) . 
This  method  was  the  four  pole  parameter  technique.  The  four  pole  param- 
eter method  had  been  widely  used  to  analyze  mechanical  systems;  therefore, 
it  was  selected  as  the  technique  to  provide  the  mathematical  verification 
of  the  beam-mounting  experiment . 

2.2  The  Four  Pole  Parameter  Technique 
2.2.1   General 

There  are  many  systems  in  nature  which  have  two  basic  quantities 
as  input  and  two  basic  quantities  as  output .  An  electric  circuit  element 
is  one  of  these  systems.  In  this  case  voltage  and  current  are  both  input 


(6) 
and  output.  All  these  systems  are  called  "four  poles"  since  four  quantities 
become  inter-related  by  some  action  within  the  system. 

Electrical  engineers  have  devised  a  technique  of  analysis  for 
circuit  elements  called  the  four  pole  parameter  technique.   Each  circuit 
element  is  considered  independently,  and  complex  circuitry  is  developed 
by  connecting  the  various  elements .  Recently  the  four  pole  parameter 
technique  has  been  applied  to  mechanical  vibration  problems  where  the 
mechanical  element  considered  has  a  force  and  velocity  as  both  input  and 
output . 

A  typical  mechanical  system  can  be  represented  in  block  diagram 
form  as  illustrated  in  Figure  II . 

Figure  II 
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F  and  V  are  the  input  force  and  velocity,  while  F  and  V  are  the 

output  force  and  velocity.   This  system  can  be  described  in  general  by 

the  following  equations; 

F  =  e   F  +  e        V  (l) 

1   £11  2    12  2  V 

V  =  e   F  +  e   V  (2) 

1    21  2    22  2  v 

e  ,  e  ,  e  ..,  and  e   are  the  four  pole  parameters  of  the  mechanical 

system. 


(7) 


Using  matrix  notation,  equations  (l)  and  (2)  become: 
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The  determination  of  the  four  pole  parameters  comprises  the  majon-portion 
of  the  technique.  Various  rules  must  also  be  developed  to  permit  combina- 
tion of  the  matrices  to  represent  mechanical  elements  in  various  multiple 
configurations . 
2.2.2   Development  of  the  Four  Pole  Parameters 

Basically  the  solution  for  the  four  pole  parameters  of  the  mechanical 
system  can  be  reduced  to  the  following  procedure  : 

1.  Investigate  the  system  and  obtain  its  performance  equations. 

2.  Solve  the  equations  imposing  the  proper  boundary  conditions. 
If  non -linearities  exist  in  the  system  it  may  be  necessary  to 
apply  linearization  techniques  at  this  step  of  the  procedure. 

3«  Reduce  the  solutions  to  the  form  of  equations  (l)  and  (2).  The 
four  pole  parameters  can  then  be  written  by  direct  comparisonu 
This  procedure  is  applied  to  several  basic  mechanical  elements  in 
Appendix  A-2. 

The  above  procedure  can  be  specifically  applied  to  the  mechanical 
systems  found  in  the  analysis  of  vibration  transmission.  Equations  (l) 
and  (2)  can  be  solved  for  the  four  pole  parameters  by  imposing  the  condition 
that  the  output,  which  is  station  2  in  Figure  II,  is  either  blocked  or 
free.  The  blocked  condition  implied  the  velocity  at  station  2  is  zero; 
the  free  condition  implies  the  force  at  station  2  is  zero. 


Fl 
€.,..  =  — —    with  station  2  blocked 
11   F2 
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Fl 
e1      -  ==—    with  station  2  free  (5) 

d       V2 

\ 

eon  =  ^r~        with  station  2  blocked  (6) 

21   F2 

e   =  —    with,  station  2  free  (7) 

The  following  definitions  can  be  applied  to  equations  (k)   through 

(7)  to  simplify  the  notation. 

F 
b  ? 

(T^)  =  Force  Transmission  =  —    with  station  2  blocked 

1 

f  V2 

(Tv)   =  Velocity  Transmission  =  —    with  station  2  free 

1 

f  V? 

(M^)  =  Transfer  Mobility  m  =£■    with  station  2  free 

b  Vl 

(M )  =  Transfer  Mobility  =  =■=■    with  station  2  blocked 
21  F2 

Equations  (k)   through  (7)  become: 
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Equation  (3)  can  now  be  written: 
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This  development  provides  a  convenient  notation  by  which  any 
system  may  be  analyzed. 
2.2-3   Combination  of  Four  Pole  Matrices 

Multiple  mechanical  elements  for  which  the  four  pole  matrices 
have  been  derived  can  be  connected  in  series .  This  is  illustrated  in 
Figure  III . 

Figure  III 


If  this  is  the  case,  the  overall  system  of  n  subsystems  S,  through  S 

can  be  represented  by  a  single  four  pole  matrix  [S  ]. 
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The  single  four  pole  matrix  [S  ]  is  defined: 

s 


[Sfl]   =    [S1]   x   [Sg]   x   [S    ]   x 


[S   -  1]  x  [Si 
n         n 


(9) 


Multiple  elements  can  also  be  connected  in  parallel  as  illustrated 


in  Figure  IV. 
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Figure  IV 


These  n  parallel  subsystems  can  be  represented  by  a  single  four  pole 

matrix  [S  ]. 
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The  derivation  of  the  matrix  [S  ]  is  more  complicated  for  the  parallel 

than  the  series  case.   The  solution  for  the  four  poles  of  [S  ]  is  included 

in  Appendix  A-3- 

2.2.4   Application  of  Four  Pole  Parameter  Technique  to  Preliminary  Experiment 

From  the  preliminary  experiment  mentioned  in  section  2.1  and  fully 
described  in  Appendix  A-l,  experimental  results  were  available  to  test  the 
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validity  of  the  four  pole  parameter  technique.  Using  this  technique  the 
isolation,  which  is  the  ratio  of  input  force  to  output  force  for  the  system, 
was  theoretically  obtained  as  a  function  of  driving  frequency.  The  experi- 
mental results  closely  agreed  with  the  values  of  isolation  determined  by  the 
four  pole  parameter  technique.  This  demonstrated  the  effectiveness  of  the 
technique  * 

2.3  The  General  System  and  Analytic  Study 
2-3-1   Description  of  the  General  System 

As  mentioned  in  Chapter  I,  the  overall  system  which  is  the  subject 
of  this  thesis  consists  of  the  mechnaical  elements  in  the  structural  path 
of  vibration  transmission  from  the  machine  source  to  the  water.   The 
block  diagram  of  this  system  is  Figure  V. 

Figure  V 
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2.3-2   The  Analytic  Study 

The  initial  step  in  any  analytic  study  is  the  enumeration  of  the 
simplifying  assumptions.   In  this  thesis  these  assumptions  are: 

1.  The  transmission  paths  in  Figure  I  can  be  considered  separately. 

2.  The  principal  contributors  to  structural  vibration  transmission 
are  the  force  and  velocity  normal  to  the  hull.  All  others  are  of 
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lesser  importance.   (This  assumption  has  been  corroborated  by  other 

investigators  (17)-) 

3.  The  mechanical  elements  are  elastic  and  linear. 

h .     The  vibrations  in  the  system  are  periodic . 
Other  assumptions  have  been  made  when  analyzing  the  individual  elements. 
These  will  be  discussed  in  deriving  the  four  poles  of  each  mechanical 
element . 

It  is  recognized  that  the  above  assumptions  will  not  permit  an 
exact  solution  of  the  vibration  transmission  characteristics  of  an  actual 
machinery  installation.  As  stated  in  Chapter  I,  however,  the  purpose  of 
this  thesis  is  not  to  accurately  predict  vibration  transmission  for  any 
particular  mechanical  system.  Rather  it  is  to  provide  a  systematic 
study  of  a  hypothetical  installation  such  that  trends  caused  by  parameter 
variations  can  be  predicted.  A  knowledge  of  these  trends  could  indicate 
design  procedures  to  reduce  vibration  transmission. 

With  the  use  of  the  simplifying  assumptions  and  the  equations  of 
dynamic  analysis,  the  four  pole  parameter  technique  can  be  applied  to 
the  general  system  described  in  Section  2.3-1*   The  dynamic  analysis 
relates  input  force  and  velocity  to  output  force  and  velocity  for  each 
element.   The  principal  task  is  the  determination  of  the  four  poles  which 
appear  in  these  equations. 

2.k     The  Four  Pole  Parameters  of  the  General  System 
2.4.1   The  Machine 

The  problem  of  mathematically  representing  the  vibration  charac- 
teristics of  an  operating  machine  is  a  very  difficult  task.  A  study  to 
evaluate  the  internal  impedance  of  a  typical  machine  has  been  made  by  the 
A.  D.  Little  Co.  (12).   In  this  report  an  ideal  machine  is  considered  to 
possess  many  resonances  and  a  complex  response  spectrum.   Even  though  in 
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some  cases  the  exciting  internal  forces  are  random,  many  machines  rotate 
and  thus  exhibit  a  form  of  periodicity  in  their  force  generation.   For 
this  reason  a  harmonic  force  generator  is  considered  as  the  vibration 
source  in  the  A.  D.  Little  report.   In  their  mathematical  model  a  harmonic 
force  generator  drives  a  multi-mass-spring  system  in  which  the  relative 
spring  stiffnesses  and  mass  sizes  determine  the  resonances.   Some  investi- 
gators have  concluded  that  this  assumed  lumped  parameter  behavior  is  valid 
at  frequencies  of  only  150  cps  or  less  (17). 

Recognizing  the  limitations  of  the  lumped  parameter  treatment, 
this  thesis  postulates  that  the  machine  can  be  represented  by  a  simple 
beam  driven  by  a  harmonic  force  generator.   A  simple  beam  by  virtue  of 
its  structural  properties  possesses  distributed  parameters  and  eliminates 
the  limitations  imposed  by  the  lumped  parameter  analysis . 

The  four  pole  parameters  for  the  beam  are  developed  in  Appendix 
A-4.1. 
2.4.2   The  Mounting 

Of  the  mechanical  systems  in  the  structural  transmission  path, 
the  mounting  is  the  one  best  approximated  by  theory.  Historically  the 
installation  of  rubber-type  resilient  mountings  was  the  first  remedial 
step  taken  to  reduce  vibrations.  Consequently,  more  mounting  theories 
have  been  developed  and  tested  than  theories  for  the  other  mechanical 
elements. 

The  theories  used  to  develop  the  four  pole  parameters  of  the  mounting 
in  the  general  system  are  those  proposed  by  D.  V.  Wright  and  A.  C.  Hagg 
(17,  7),  A.  0.  Sykes  (13,  Ik) ,    and  J.  C.  Snowden  (18). 

The  first  mounting  considered  in  the  general  system  is  the  standard 
Portsmouth  mount,  which  is  found  in  many  shipboard  installations  today. 
The  describing  four  poles  for  the  mount  are  derived  in  Appendix  A-4.2 
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according  to  the  theory  developed  "by  Wright  and  Hagg. 

The  second  mounting  considered  in  the  general  system  is  one 
constructed  of  a  resilient  rubber-type  material  similar  to  that  bearing 
the  trade  name  of  Thiokol  R.  D.   The  four  pole  parameters  for  this  mount 
are  developed  following  the  concepts  of  Sykes  and  Snowden  in  Appendix 
A-4.2. 
2.4-3   The  Foundation 

In  most  practical  installations,  machinery  foundation  structures 
are  combinations  of  simple  beams.   This  fact  naturally  suggests  the  use 
of  beam  vibration  analysis  to  describe  the  foundation  element.   This 
analysis  has  been  used  by  a  number  of  investigators  as  a  basic  step  in  the 
description  of  the  dynamics  of  foundation  vibration  (17>  20 ) . 

For  the  purpose  of  this  investigation  the  foundation  was  simplified 
to  be  a  beam  resting  on  two  columns.   The  describing  four  poles  for  this 
mechanical  system  are  derived  in  Appendix  A-4.3* 
2.k.k       The  Hull 

The  hull  of  a  ship  is  a  complex  structure.   Unlike  other  elements 
in  the  structural  vibration  transmission  path,  it  cannot  be  represented  by 
any  single  mechanical  element  or  simple  combination  of  elements .   Experi- 
mental results  are  necessary  to  describe  the  response  of  typical  hull 
sections  (7,  11,  17,  20). 

A  mathematical  model  devised  by  D.  V.  Wright  has  provided  a 
reasonable  approximation  of  the  vibration  characteristics  of  a  ship's  hull. 
This  model  is  used  to  obtain  the  four  poles  of  the  hull  section  in  Appendix 
A-k.k. 
2.4.5   The  Water 

The  water  is  considered  only  in  its  effect  as  an  attenuation 
element.  Any  entrained  mass  effect  is  included  as  part  of  the  hull.  The 
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water  presents  an  effective  specific  resistance  to  the  propagation  of 
logitudinal  pressure  waves.   This  specific  resistance  is  dependent  on  the 
density  of  sea  water  and  the  velocity  of  sound  in  that  medium. 

The  describing  four  poles  are  derived  in  Appendix  A-k.5. 
2-5  Solution  of  the  Problem 
2.5-1   General 

The  derivation  of  an  expression  for  vibration  transmission 
through  the  mechanical  elements  of  the  general  system  is  the  key  to 
achieving  the  purpose  of  this  thesis.  The  solution  of  this  problem 
involves  the  combination  of  the  five  basic  subsystems  described  in  Section 
2.4  into  a  single  complex  system.   The  four  pole  parameters  of  the  combined 
system  are  then  obtained  according  to  the  equations  of  combination  explained 
in  Section  2.2.3-   Ike  vibration  transmission  is  defined  as  a  force  ratio. 
The  force  transmitted  into  the  water  is  considered  the  significant  noise 
producing  force  in  the  system.   Therefore,,  the  vibration  transmission  is 
defined  as  the  ratio  of  the  force  into  the  water  to  the  input  force . 

Since  many  of  the  four  poles  in  the  system  are  functions  of  fre- 
quency, the  final  ratio  will  be  a  function  of  frequency.   Evaluation  of  ';  . 
this  ratio  at  sufficient  frequencies  will  give  the  desired  transmission 
versus  frequency  data.   To  conform  to  convention  transmission  is  converted 
to  decibels . 

The  combination  of  the  general  system  fourjppoles  and  a  sample 
calculation  at  one  value  of  frequency  are  included  in  Appendix  A-6. 
2.5.2   Computer  Application 

This  thesis  could  not  have  been  completed  without  the  computational 
assistance  of  the  IBM  709  digital  computer.   The  solution  for  vibration 
transmission  at  each  frequency  required  a  calculation  of  length  comparable 
to  that  of  Appendix  A-6.   This  study  required  thousands  of  such  calcula- 
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tions.  The  matrix  notation  of  the  four  pole  parameter  technique  was 
easily  adapted  to  digital  computer  methods . 

A  computer  program  was  written  to  solve  the  combination  of  five 
four  pole  matrices.  Provisions  were  made  to  change  any  of  the  mechanical 
elements  in  the  system  of  Figure  V.  A  complete  sequence  of  this  computer 
program  logic  is  explained  in  Appendix  A~7«  A  most  convenient  feature  of 
the  program  was  the  use  of  the  cathode  ray  tube  oscilloscope  to  graphically 
present  the  results.   The  scope  photograph  eliminated  the  necessity  of 
manual  plotting. 
2.6  Systematic  Study  to  Determine  Trends 

With  the  selection  of  the  general  system  of  Figure  V  and  the 
mechanical  elements  in  Section  2.k   to  represent  this  system,  the  problem  to 
be  considered  was  established.   The  four  pole  parameter  technique  described 
in  Section  2.2  provided  an  analytical  means  to  study  the  problem.   The 
final  phase  of  the  study,  however,  was  the  most  important.  It  was  to 
devise  a  consistent,  systematic  method  of  changing  the  general  system 
such  that  many  different  cases  might  be  evaluated  and  compared.   This 
comparison  and  evaluation  formed  the  basis  for  determining  the  results, 
conclusions,  and  recommendations  of  this  thesis. 

A  detailed  description  of  the  systematic  changes  to  the  general 
system  is  given  in  Appendix  A~5* 
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CHAPTER  III 
RESULTS 

3.1  The  Photo  Plots 

The  results  of  the  analytic  study  of  vibration  transmission  in 
shipboard  installations  are  presented  in  the  form  of  forty-one  photo  plots. 
These  transmission  versus  frequency  plots  are  included  in  Appendix  B  with 
one  plot  for  each  of  the  cases  mentioned  in  Appendix  A-5'7* 

3.2  Statements  of  Results 

The  following  statements  of  results  are  deduced  from  a  study  of 
the  photo  plots.   To  preserve  the  systematic  procedure  of  the  study,  each 
case  is  examined  separately.  For  detailed  descriptions  of  the  systems 
studied,  refer  to  Appendix  A-5 • 

Case  1.   This  was  the  basic  case  involving  System  1.   It  was 

selected  as  the  basis  for  comparison  since  the  elements  comprising 

System  1  were  typical  of  those  presently  found  in  shipboard 

installations . 

Case  2 .   The  System  was  the  same  as  Case  1,  except  for  the 

following : 

a.  Mount  damping  increased. 

The  first  transmission  peak  was  lowered. 

Case  3«   The  System  was  the  same  as  Case  1,  except  for  the  following: 
a.  Mount  stiffness  increased. 

The  transmission  over  the  frequency  range  studied  was  increased. 

Case  k.      The  System  was  the  same  as  Case  1,  except  for  the  following: 
a.  Mount  stiffness  decreased. 

The  transmission  over  the  frequency  range  studied  was  lowered  by 

approximately  10  db. 

Case  5 •  The  System  was  the  same  as  Case  1,  except  for  the 
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following : 

a.  Thiokol  R.  D.  mount  replaced  Portsmouth  mount. 
The  first  transmission  peak  was  lowered  while  the  second  peak 
was  removed. 
Case  6.  The  System  was  the  same  as  Case  1,  except  for  the  following: 

a.  Hull  damping  was  doubled. 
There  was  no  change  from  Case  1. 
Case  7»  The  System  was  the  same  as  Case  1,  except  for  the  following : 

a.  Mount  damping  increased. 

b.  Mount  stiffness  decreased. 
The  results  were  the  same  as  Case  k. 

Case  8.  The  System  was  the  same  as  Case  1,  except  for  the  following: 

a.  Thiokol  R.  D.  mount  replaced  Portsmouth  mount. 

b.  Hull  damping  was  doubled. 
The  results  were  the  same  as  Case  5. 

Case  9»  1he  System  was  the  same  as  Case  1,  except  for  the  following: 

a.  Mount  damping  increased. 

b.  Mount  stiffness  decreased. 

c .  Hull  damping  increased  twenty  times . 

The  transmission  was  two  or  three  decibels  below  that  of  Case  7  a"t 
the  lower  frequencies. 

Case  10.   The  System  was  the  same  as  Case  1,  except  for  the 
following : 

a.  Viscous  foundation  damping  was  added. 
The  initial  transmission  slope  was  changed.  The  transmission  was 
lowered  at  frequencies  greater  than  1700  cps. 
Case  11.   The  System  was  the  same  as  Case!  1,  except  for  the 
following : 
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a.  Thiokol  R.  D.  and  Portsmouth  mounts  in  parallel. 
The  transmission  was  increased  at  frequencies  greater  than  1000  cps. 
Case  12 ♦  The  System  was  the  same  as  Case  1,  except  for  the 
following : 

a.  Water  resistance  varied  with  frequency. 
The  transmission  was  lowered  at  frequencies  greater  than  1700  cps. 
Case  13 ♦  The  System  was  the  same  as  Case  1,  except  for  the 
following : 

a.  Viscous  foundation  damping  was  added. 

b.  Mount  stiffness  decreased. 

The  first  transmission  peak  was  shifted  from  21  cps  to  11  cps. 
Transmission  was  decreased  over  the  entire  frequency  range  "by 
at  least  10  db. 

Case  l4.  The  System  was  the  same  as  Case  1,  except  for  the 
following : 

a.  Thiokol  R.  D.  foundation  damping  was  added. 
Transmission  was  lowered  appreciably  at  frequencies  above  800  cps. 
Case  15.   The  System  was  the  same  as  Case  1,  except  for  the 
following  2 

a.  Thiokol  R.  D.  and  Portsmouth  mounts  in  parallel. 

b.  Portsmouth  mount  stiffness  decreased. 

The  results  were  the  same  as  Case  k   to  500  cps.  Above  500  cps, 
the  results  were  the  same  as  Case  11. 

Case  16.   The  System  was  the  same  as  Case  1,  e  cept  for  the 
following : 

a.  Hull  damping  increased,  maximum. 
The  transmission  was  reduced  over  the  entire  frequency  range. 
Case  17.  The  System  was  the  same  as  Case  1,  except  for  the 
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following : 

a.  Thiokol  R.  D.  mount  replaced  Portsmouth  mount. 

b.  Thiokol  R.  D.  foundation  damping  was  added. 

The  first  and  second  transmission  peaks  were  removed.  The 
transmission  was  reduced  at  frequencies  above  800  cps. 
Case  18.  The  System  was  the  same  as  Case  1,  except  for  the 
following : 

a.  Thiokol  R.  D.  foundation  damping  was  added. 

b.  Mount  stiffness  decreased. 

The  vibration  transmission  was  reduced  10  db.  in  the  10  to  800  cps 
range.  Above  800  cps,  the  transmission  was  further  decreased. 
Case  19*  The  System  was  the  same  as  Case  1,  except  for  the 
following : 

a.  Thiokol  R.  D.  foundation  damping  was  added. 

b.  Mount  stiffness  decreased. 

c.  Mount  damping  increased. 

The  results  were  the  same  as  Case  18,  except  the  first  peak  was 

removed. 

Case  20.  This  was  the  basic  case  involving  System  2.  Foundation 

mass  increase  was  the  only  change  from  System  1.  The  transmission 

was  lowered  at  frequencies  from  10  to  1000  cps.  The  second  and 

third  transmission  peaks  were  split. 

Case  21.  The  System  was  the  same  as  Case  20,  except  for  the 

following : 

a.  Mount  stiffness  decreased. 
The  transmission  over  the  frequency  range  studied  was  lowered  by 
approximately  10  db. 
Case  22 »  The  System  was  the  same  as  Case  20,  except  for  the 
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following  i 

a.  Thiokol  R.  D.  mount  replaced  Portsmouth  mount. 
The  first  transmission  peak  was  lowered;  the  second  split  peak  was 
removed;  and  half  of  the  third  split  peak  was  removed. 
Case  23 ■  The  System  was  the  same  as  Case  20,  except  for  the 
following : 

a.  Mount  damping  increased. 

b.  Mount  stiffness  decreased. 

c .  Hull  damping  increased  twenty  time  s . 

The  first  transmission  peak  was  shifted  below  10  cps,  and 
transmission  over  the  entire  frequency  range  was  lowered  by  10  db. 
Case  2h.      The  System  was  the  same  as  Case  20,  except  for  the 
following : 

a.  Viscous  foundation  damping  was  added. 
The  split  transmission  peak  at  3^-0  cps  was  lowered,  and  the 
transmission  from  1000  to  10,000  cps  was  lowered  by  5  to  10  db. 
Case  25 .  The  System  was  the  same  as  Case  20,  except  for  the 
following : 

a.  Thiokol  R.  D.  mount  replaced  Portsmouth  mount. 

b.  Thiokol  R.  D.  foundation  damping  was  added. 

The  first  transmission  peak  was  removed  while  the  second  and  third 
peaks  were  lowered.   The  transmission  was  considerably  reduced  at 
frequencies  above  300  cps . 

Case  26.   The  System  was  the  same  as  Case  20,  except  for  the 
following : 

a.  Thiokol  R.  D.  foundation  damping  was  added. 

b.  Mount  stiffness  decreased. 

The  first  transmission  peak  was  shifted  below  10  cps.  With  the 
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exception  of  the  second  peak,  the  vibration  transmission  was 
considerably  reduced  over  the  entire  frequency  range. 
Case  30»  This  was  the  basic  case  involving  System  3*  Foundation 
mass  increase  was  the  only  change  from  System  2.  The  transmission 
was  considerably  lowered  over  the  entire  range  except  for  peaks  at 
800  and  ^000  cps. 

Case  31 ♦  The  System  was  the  same  as  Case  30,  except  for  the 
following : 

a.  Mount  stiffness  decreased. 
The  first  transmission  peak  was  removed,  and  the  transmission 
reduced  by  5  to  10  db.  over  the  frequency  range. 
Case  32.  The  System  was  the  same  as  Case  30,    except  for  the 
following  j 

a.  Thiokol  R.  D.  mount  replaced  Portsmouth  mount. 
The  first,  second,  and  third  transmission  peaks  were  removed. 
Case  33 ■  The  System  was  the  same  as  Case  30,  except  for  the 
following : 

a.  Mount  damping  increased. 

b.  Mount  stiffness  decreased. 

c.  Hull  damping  increased  twenty  times. 
The  results  were  the  same  as  Case  31» 

Case  3^-»  The  System  was  the  same  as  Case  30,  except  for  the 
following : 

a.  Viscous  foundation  damping  was  added. 
The  first  transmission  peak  was  split;  the  peaks  at  800  and  2500 
cps  were  lowered.  The  peak  at  4000  cps  was  removed. 
Case  35 ■   The  System  was  the  same  as  Case  30,  except  for  the 
following : 
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a.  Thiokol  R,  D.  mount  replaced  Portsmouth  mount. 

b.  Thiokol  R.  D.  foundation  damping  was  added. 

With  the  exception  of  the  transmission  peak  at  200  cps,  the 

vibration  transmission  was  considerably  lowered  over  the  range  of 

frequencies. 

Case  36.  The  System  was  the  same  as  Case  30,  except  for  the 

following : 

a.  Thiokol  R.  D.  foundation  damping  was  added. 

b.  Mount  stiffness  decreased. 

The  first  transmission  peak  was  removed,  and  the  transmission 

reduced  by  5  to  10  db.  up  to  500  cps.  Above  500  cps,  the 

transmission  was  further  reduced. 

Case  50.  This  was  the  basic  case  involving  System  h.     The 

machine  and  foundation  masses  were  ten  times  the  values  in 

System  1.  This  System  was  used  as  the  basis  of  comparison  for 

larger  machinery  installations. 

Case  51 •  The  System  was  the  same  as  Case  50,  except  for  the 

following : 

a.  Mount  stiffness  decreased. 
The  first  transmission  peak  was  shifted  and  lowered.  The 
transmission  was  reduced  by  5  to  10  db.  over  the 
Case  52 •  The  System  was  the  same  as  Case  50,  except  for  the 
following : 

a.  Thiokol  R.  D.  mount  replaced  Portsmouth  mount. 
The  first  transmission  peak  was  removed,  and  the  transmission  was 
lowered  by  5  to  10  db.  over  the  entire  range. 
Case  53-  The  System  was  the  same  as  Case  50,  except  for  the 
following : 
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a.  Mount  damping  increased. 

b.  Mount  stiffness  decreased. 

c.  Hull  damping  increased  twenty  times. 

The  first  transmission  peak  was  considerably  lowered.  The 
transmission  was  lowered  over  the  entire  range  by  approximately 
10  db. 

Case  5^«  1^e  System  was  the  same  as  Case  50,  except  for  the 
following : 

a.  Viscous  foundation  damping  was  added. 
The  first  transmission  peak  was  lowered;  the  third  peak  was 
eliminated.  The  transmission  was  lowered  by  5  "to  10  db.  from  1000 
to  10,000  cps. 

Case  55 •  The  System  was  the  same  as  Case  50,  except  for  the 
following : 

a.  Thiokol  R.  D.  mount  replaced  Portsmouth  mount. 

b.  Thiokol  R.  D.  foundation  damping  was  added. 

The  first  transmission  peak  was  removed,  and  the  transmission  was 
reduced  by  10  to  15  db.  to  300  cps.  Above  300  cps,  it  was  reduced 
by  20  to  kO   db. 

Case  56.      The  System  was  the  same  as  Case  50,  except  for  the 
following : 

a.  Thiokol  R.  D.  foundation  damping  was  added. 

b.  Mount  stiffness  decreased. 

The  first  transmission  peak  was  lowered  and  shifted  to  12  cps. 
The  second  peak  was  reduced  by  35  db.  At  higher  frequencies  the 
transmission  was  extremely  low. 

Case  57*  The  System  was  the  same  as  Case  50,  except  for  the 
following  1 
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a.  Thiokol  R.  D.  foundation  damping  was  added. 

b.  Mount  stiffness  decreased. 

c.  Mount  damping  increased. 

The  first  transmission  peak  was  removed.  The  second  peak  was 
reduced  by  35  clb.  At  higher  frequencies  the  transmission  was 
extremely  low. 
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CHAPTER  IV 
DISCUSSION  OF  RESULTS 
k.l     Preliminary  Discussion 

The  most  important  single  result  of  this  thesis  is  the  technique:- 
of  analytic  description  of  a  complex  shipboard  machinery  installation. 
This  technique  formulates  a  workable  analytic  tool  for  the  study  of 
vibration  transmission  characteristics  of  these  complex  systems.  In 
essence  this  technique  is  the  four  pole  parameter  method.  It  is  recom- 
mended that  this  method  of;  analysis  be  considered  in  future  design  inves- 
tigations . 

k.2     Systematic  Changes 
4.2.1   Mounting 

The  properties  of  the  resilient  mounting  had  a  profound  effect 
on  the  vibration  transmission.  A  study  of  the  photo  plots  showed  that 
stiffness  in  the  Portsmouth  mounting  had  a  dominant  iffect  in  the  cases 
in  which  it  was  altered.  When  stiffness  was  increased,  the  transmission 
at  all  frequencies  was  increased;  when  stiffness  was  decreased,,  the 
transmission  at  all  frequencies  was  lowered  by  10  db.  It  is  recommended 
that  the  design  engineer  make  stiffness  as  low  as  possible  and  yet 
consistent  with  shock,  alignment,  and  other  factors  influencing  mount 
design. 

Damping  in  the  Portsmouth  mounting  was  not  nearly  as  influential 
as  mount  stiffness.  However,  increased  damping  was  effective  in  decreasing 
the  transmission  peaks  at  lower  frequencies.  It  is  recommended  that  an 
effort  be  made  to  increase*  the  damping  of  the  Portsmouth  mount  by  increas- 
ing the  material  damping  or  by  Introducing  an  external  damping  device  in 
parallel  with  the  mount. 

The  replacement  of  the  Portsmouth  mounting  by  a  synthetic  rubber 
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mounting  of  Thiokol  R.  D.  material  showed  interesting  results.  The 
Thiokol  R.  D.  mount  lowered  the  first  transmission  peak  and  eliminated  the 
second.  Therefore,  the  use  of  the  Thiokol  R.  D.  mount  appears  desirable 
in  systems  with  high  noise  transmission  at  low  frequencies  or  in  systems 
in  which  low  frequency  noise  suppression  is  paramount.  The  stiffness  of 
Thiokol  R.  D.  increases  exponentially  with  frequency  (18).  This  decreases 
the  effectiveness  of  Thiokol  R.  D.  in  suppressing  vibration  transmission 
at  higher  frequencies.  The  development  of  a  material  to  correct  this 
undesirable  trend  of  Thiokol  R.  D.  is  recommended. 
4.2.2   Foundation 

Of  the  various  foundation  parameters  the  mass  most  affected  the 
vibration  transmission.   Three  different  foundation  masses  were  investi- 
gated with  the  result  that  vibration  transmission  was  lowered  as  the 
foundation  mass  increased.  Within  the  limitations  of  weight  and  space, 
foundations  of  vibrating  machinery  should  be  made  as  massive  as  possible. 

Viscous  foundation  damping  and  damping  utilizing  the  Thiokol  R.  D. 
material  were  studied.  Foundation  damping  was  effective  at  frequencies 
above  1000  cps  for  the  vispous  type  and  above  800  cps  for  the  Thiokol  R.  D. 
type.  However,  the  vibration  transmission  was  reduced  more  by  the  Thiokol 
R.  D.  type. 

Foundation  damping  is  recommended  in  installations  where  high 
frequency  noise  transmission  is  a  problem. 
4.2-3   Hull 

The  effects  of  increasing  the  hull  damping  first  became  apparent 
in  the  mathematical  model  used  when  the  damping  was  increased  by  an  order 
of  twenty.  As  hull  damping  was  further  increased,  the  vibration  transmission 
was  reduced  at  all  frequencies. 

These  results  are  at  best  qualitative.  In  this  model  lumped 


(28) 
parameter  damping  was  used,  and  no  attempt  was  made  to  study  the  effects 
of  distributed  damping.  In  practice  the  latter  would  be  the  type  applied 
to  the  hull.  A  study  to  quantitatively  evaluate  the  effects  of  distributed 
hull  damping  would  be  useful  In  the  extension  of  an  analytic  study. 
4.2 A   Combined  Changes 

In  a  series  of  cases  studied  the  Portsmouth  mount  stiffness  was 
decreased,  the  mount  damping  increased,  and  hjoll  damping  Increased  twenty 
times.  This  was  done  with  the  four  systems  studied.  In  general,  as  the 
masses  in  the  systems  were  increased,  the  vibration  transmission  was 
reduced.  The  results  of  this  series  study  were  compared  with  those  cases 
In  which  mount  stiffness  alone  was  changed.  This  comparison  showed  that 
the  change;  in  mount  stiffness  overrode  the  changes  in  mount  damping  and 
hull  damping. 

A  single  case  of  decreased  Portsmouth  mount  stiffness  and  increased 
mount  damping  showed  that  the  vibration  transmission  was  lowered  by  exactly 
the  same  amount  as  the  case  in  which  stiffness  alone  was  decreased.  This 
further  supported  the  conclusion  that  mount  stiffness  was  one  of  the  most 
important  factors  to  be  considered  in  any  vibration  transmission  .study. 

When  a  mount  of  Thiokol  R.  D.  was  placed  in  parallel  with  a 
Portsmouth  mount,  the  transmission  reduction  properties  were  not  improved. 
The  addition  of  the  Thiokol  R.  D.  mount  actually  increased  the  transmission 
at  frequencies  greater  than  1000  cps.  This  case  was  again  studied  with 
the  Portsmouth  mount  stiffness  decreased.  Results  showed  that  this  case 
was  not  as  good  as  that  in  which  the  single  Portsmouth  mount  with  decreased 
stiffness  was  examined.  Therefore,  it  is  concluded  that  placing  the 
Thiokol  R.  D.  mounting  in  parallel  with  a  Portsmouth  mounting  would  not  be 
advantageous  in  reducing  vibration  transmission.  If  parallel  mountings 
are  to  be  used,  as  recommended  by  Snowden  in  Reference  (18),  care  must  be 
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used  In  the  selection  of  the  high  damping  synthetic  rubber  material. 

In  another  series  of  cases  the  Portsmouth  mount  was  replaced  by 
a  Thiokol  R.  D.  mount  and  Thiokol  R.  D.  damping  was  applied  to  the  founda- 
tion. This  was  done  with  the  four  systems  studied.  As  the  masses  in  the 
system  were  increased,  the  vibration  transmission  decreased  markedly.  In 
any  one  system  the  combination  of  Thiokol  R.  D.  mount  and  Thiokol  R.  D. 
foundation  damping  was  most  effective  in  reducing  transmission.  It  Is 
recommended  that  this  combination  be  investigated  in  an  experimental 
installation. 

A  Portsmouth  mount  with  decreased  stiffness  and  Thiokol  R.  D. 
foundation  damping  comprised  the  changes  examined  in  another  series  of 
cases.  Again  this  was  done  with  the  four  systems.  The  vibration 
transmission  was  significantly  decreased  as  the  masses  in  the  system  were 
increased.  Within  a  particular  system  this  combination  was  less  effective 
at  frequencies  of  less  than  kO   cps  than  the  Thiokol  R.  D.  mount  and  Thiokol 
R.  D.  foundation  damping  combination.  However,  It  was  superior  at 
frequencies  higher  than  ^0  cps.  This  combination  should  also  be  investi- 
gated in  an  experimental  installation. 

A  Portsmouth  mount  with  decreased  stiffness  and  increased  damping 
and  Thiokol  R.  D.  foundation  damping  comprised  the  changes  in  the  final 
series  examined.  The  only  transmission  change  .  from  the  previous  series 
was  the  elimination  of  the  first  transmission  peak. 
4.3  General  Considerations 

The  results  of  the  systematic  study  largely  verified  existing 
concepts  in  reducing  vibration  transmission  in  shipboard  machinery 
installations.  This  would  indicate  that  the  assumptions  and  the  four  pole 
parameter  method  of  analysis  were  sufficient  to  establish  the  vibration 
trends  and  characteristics  proposed  at  the  beginning  of  this  thesis. 
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As  new  ideas  for  reducing  vibration  transmission  are  introduced,  it  is 
proposed  that  they  be  investigated  by  means  of  the  technique  of  this 
thesis. 

The  four  pole  parameter  technique  enables  use  of  empirical  data 
obtained  independently  for  each  mechanical  element  in  the  vibration  path. 
Empirical  data  can  be  used  for  elements  which  defy  analytic  description. 
It  would  appear  that  the  use  of  empirical  data  for  the  element  four  poles 
would  more  closely  approximate  the  vibration  characteristics  of  an  actual 
system. 

An  interesting  extension  of  this  thesis  would  be  an  optimization 
study  to  determine  the  best  combination  of  elements  in  the  vibration 
transmission  paths.  Since  shock  specifications  also  govern  the  element 
selection,  the  optimization  must  include  shpck  as  well  as  vibration 
reduction  considerations. 

The  construction  of  a  model  of  the  approximate  system,  as 
represented  in  Figure  A-3  of  Appendix  A,  would  permit  an  experimental 
verification  of  the  four  pole  parameter  technique. 
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CHAPTER  V 
CONCLUSIONS 
As  logically  developed  in  Chapter  IV  the  following  statements 
are  listed  as  the  conclusions  of  this  thesis : 

1.  The  four  pole  parameter  technique  provides  a  sufficiently 
powerful  method  of  analyzing  vibration  transmission.  Furthermore, 
the  four  pole  parameter  technique  is  easily  adapted  to  a  digital 
computer  program. 

2.  For  the  purpose  of  vibration  analysis  a  complex  machinery 
installation  can  be  represented  by  a  relatively  simplified  system. 
This  representation  is:  sufficient  to  establish  vibration  trends 
and  characteristics. 

3«  The  stiffness  of  the  resilient  mounting  has  a  dominant  effect 
on  the  vibration  transmission  at  all  frequencies. 
k.     The  only  effect  of  increased  damping  in  the  Portsmouth  mount 
was  to  lower  the  first  vibration  transmission  peak. 

5 .  A  mount  fabricated  from  a  high  damping  synthetic  rubber  such 
as  Thiokol  R.  D.  is  highly  effective  in  reducing  vibration 
transmission  at  low  frequencies.  The  same  mount  is  not  as  effective 
as  the  Portsmouth  mount  at  higher  frequencies. 

6.  As  the  mass  of  the  foundation  is  increased,  the  vibration 
transmission  is  lowered. 

7.  Foundation  damping  was  effective  at  frequencies  above  1000 
cps  for  the  viscous  type  and  800  cps  for  the  Thiokol  R.  D.  type. 
The  Thiokol  R.  D.  type  was  more  effective. 

8.  Increasing  hull  damping  decreases  vibration  transmission  at '.all 
frequencies . 

9»  Placing  a  Thiokol  R.  D.  type  mounting  tin  parallel  with  a 
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Portsmouth  mount  as  suggested  in  Reference  (l8)  actually  bad  a 
detrimental  effect  on  vibration  reduction. 

10.  Replacing  a  normal  mount  by  a  Thiokol  R.  D.  mount  and  applying 
Thiokol  R.  D.  damping  to  the  foundation  were  most  effective  in 
reducing  vibration  transmission. 

11.  Decreasing  the  stiffness  of  a  Portsmouth  mount  and  applying 
Thiokol  R.  D.  damping  to  the  foundation  were  superior  to  the 
system  mentioned  in  the  previous  conclusion  in  reducing  vibrations 
at  frequencies  above  ^0  cps# 

12.  A  Portsmouth  mount  with  decreased  stiffness  and  increased 
damping  and  Thiokol  R.  D.  damping  applied  to  the  foundation  were 
the  most  effective  combination  in  overall  transmission  reduction. 
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CHAPTER  VI 
RECOMMENDATIONS 
As  stated  or  inferred  in  Chapter  IV  the  following  are  the 
recommendations  of  this  thesis: 

1.  The  four  pole  parameter  method  of  analysis  should  be  considered 
in  future  design  investigations  involving  vibration  transmission. 

2.  As  new  ideas  for  reducing  vibration  transmission  are  introduced, 
they  should  be  investigated  by  means  of  the  technique  of  this  thesis 

3.  Empirical  data  should  be  used  for  element  four  pole  parameters 
to  more  closely  approximate  the  vibration  characteristics  of  an 
actual  system. 

k»     An  optimization  study  to  determine  the  best  combination  of 
elements  in  the  vibration  transmission  paths  should  be  conducted. 
It  should  include  shock  as  well  as  vibration  considerations. 
5 •  A  model  of  the  approximate  system  of  this  thesis  should  be 
constructed  to  permit  experimental  verification  of  the  four  pole 
parameter  technique . 

6.  Mount  stiffness  should  be  as  low  as  possible  and  yet  consis- 
tent with  shock,  alignment,  and  other  factors  influencing  the 
design. 

7.  The  material  damping  of  the  Portsmouth  mount  should  be 
increased,  or  an  external  damping  device  placed  in  parallel  with 
the  mount. 

8.  The  Thiokol  R.  D.  mount  should  be  used  in  systems  with  high 
noise  transmission  at  low  frequencies  or  in  systems  in  which  low 
frequency  noise  suppression  is  paramount. 

9.  A  high  damping  mounting  material  should  be  developed  in 
which  stiffness  remains  constant  or  decreases  with  frequency. 
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10.  Within  the  limitations  of  weight  and  space,  foundations  of 
vibrating  machinery  should  be  as  massive  as  possible. 

11.  Foundation  damping  should  be  utilized  in  installations  where 
high  frequency  noise  transmission  is  a  problem. 

12.  A  study  to  quantitatively  evaluate  the  effects  of  distributed 
hull  damping  should  be  conducted. 

13.  Caution  should  be  exercised  in  the  selection  of  a  high  damping 
synthetic  rubber  mount  to  be  placed  in  parallel  with  a  standard 
mount. 

1^-.  The  effects  of  replacing  a  normal  mount  by  a  Thiokol  R.  D. 
mount  and  applying  Thiokol  R.  D.  damping  to  the  foundation  should 
be  investigated  in  an  experimental  installation. 

15.  The  effects  of  decreasing  the  stiffness  of  a  Portsmouth  mount 
and  applying  Thiokol  R.  D.  damping  to  the  foundation  should  be 
investigated  in  an  experimental  installation. 

16.  The  effects  of  decreasing  the  stiffness  and  increasing  the 
damping  of  a  Portsmouth  mount  and  applying  Thiokol  R.  D.  damping 
to  the  foundation  should  be  investigated  in  an  experimental 
installation. 
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APPENDIX  A 


DETAILS  OF  PROCEDURE 
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A-l  The  Be am-Mount ing  System 
A-l .  1   Exp  e  r  iment 

This  experiment  was  conducted  to  determine  the  isolation  charac- 
teristics of  a  typical  Admiralty  mounting.   The  experimental  system  consisted 
of  a  fifty  pound  beam  driven  by  a  harmonic  vibration  generator.  The  beam 
was  supported  by  the  mounting  which  rested  on  a  "rigid"  foundation.  The 
components  and  instrumentation  are  illustrated  in  Figure  A-l. 

The  foundation  system  was  not,  in  fact,  rigid.  However,  its  natural 
frequency  was  much  lower  than  the  rigid  body  mode  frequency  of  the  beam-:1. 
mounting  system.  This  insured  no  resonance  interaction,  and  the  foundation 
mass,  m_  in  Figure  A-l,  responded  in  a  rigid  manner.  The  vibration 
generator  imparted  a  velocity  to  the  mass  hl.  The  velocity  V_  into  the 
mounting  of  compliance  C.  and  the  velocity  V-  at  the  free  side  of  the  beam 
were  picked  up  by  the  barium  titanate  accelerometers,  A  and  A_^  respec- 
tively. The  isolation  in  this  case  is  the  ratio  of  V  to  V_ .  Both  this 
and  the  driving  frequency  were  recorded,  and  the  isolation  was  converted 
to  decibels .  This  was  repeated  over  the  frequency  range  of  10  to  3000  cps . 

To  obtain  a  more  realistic  model  in  which  the  beam  was  driven 
from  above,  as  a  machine  normally  would  be,  use  of  the  transmissibility 
theorem  was  necessary.  This  stated  that  the  transmission  (or  isolation) 
in  one  direction  for  linear  elastic  systems  is  equivalent  to  the  force 
transmission  in  the  opposite  direction: 

V  F 

VI  F2 

This  new  system  could  be  represented: 
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The  physical  constants  of  the  system  of  the  experiment  were  as  follows 
Beam:  Length:  L  =  60  in. 

Thickness:  t  =  -75  in» 
Weight:  W_  =  50  lb. 
Frequencies : 


> ■-*  J- 

n    n  v  : 


EI 
co„  =  /i   /  — 

"B 


where:  E  =  30  x  10  psi 

I  =  moment  of  inertia.  =  -=-r-  — — 

12  pL 

p  =  mass  density 

Describing  Functions: 

<Zf,  (— )  and  ^  given  in  Reference  (l6) 
Hi  v2       n 

Mounting :  Compliance :  C  =  10rn  tt1 

A-1.2   Theory 

The  beam-mounting  system  was  represented  schematically: 
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. S  ]  was  the  matrix  of  four  pole  parameters  for  the  beam  and 


[S^]  for  the  mounting. 
Then: 
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or  by  individual  steps : 
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The  matrices  in.  terms  of  their  four  poles  were  expressed: 
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As  developed  in  Section  2.2.2,  the  general  four  pole  parameter 
matrix  is: 
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In  this  proof  of  the  experiment  the  mounting  was  treated  as  a 
simple  spring.   From  Appendix  A-2,  the  four  poles  for  a  simple  spring  are? 
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The  compliance  of  the  spring  is  C  ,  while  the  mobility  is  jwC  . 

The  beam  was  analyzed  following  the  techniques  of  Ruzicka  and 
Cavanaugh  (h) .     With  negligible  compression  assumed  in  the  transverse 
direction  and  from  the  transmissibility  theory,  the  four  poles  of  the  beam 
were  evaluated  : 
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where  M(— )  was  the  mobility  at  the  center  of  the  beam.   This  was  taken 
from  Reference  [k)i 
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The  matrix  [S- ]  for  the  beam  was  then  written: 
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Equations  (5)  and  (7)  were  substituted  in  equation  (l): 
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To  check  the  experimental  results  a  plot  of  isolation  I,  which  is 

defined  as  =— ,  versus  frequency  was  desired.   Equation  (8)  was  rewritten: 
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x  =  (JWC^)  F3  +  V3 


Station  3  in  the  experiment  was  a  rigid  support;  therefore,  V  =  0 
The  isolation  was  obtained: 
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The  mobility  of  the  beam  is  expressed  in  Equation  (6). 
A-l.3   Comparison  of  Experiment  and  Theory 

The  curves  of  isolation  in  Figure  A-2  show  that  the  four  pole 
parameter  technique  provided  excellent  theoretical  verification  of  the 
experimentally  measured  values  of  isolation.   These  results  indicated  the 
value  of  the  four  pole  analysis  and  suggested  the  use  of  it  in  the  study 
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FIGURE    A-2 
PRELIMINARY     EXPERIMENT 
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of  more  complex  systems . 

A-2  Development  of  Four  Poles 

A-2.1   Mass  Four  Poles 


The  performance  equations  are 


Vl   V2 


because  the  mass  is  a  rigid  body  and 


dV     dv2 
Fi  "  F2  =  m  dt"  =  m  dT- 


(1) 
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from  the  application  of  Newton's  Law. 

Since  the  applied  force  is  usually  considered  as  harmonic,  the 
other  force  and  the  velocities  will  also  be  harmonic  functions . 

F  =  F   e^Wt    F  =  F   e^ 
1    10  2    20  e 

1    10         V2    20 
F-n,  ^_n,  ^o'  an(^  ^on  are  time  independent  but  may  be  complex. 
Substituting  the  expression  for  V_ ,  Equation  (2)  becomes: 
F]_  =  f2  +  ( jm)   v2 

while  (l)  remains 

v1  =  v2 
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By  inspection  the  four  pole  parameters  are : 
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A-2  .2   Spring  Four  Poles 
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The  performance  equations  are 
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because  the  spring  possesses  no  inertia  and 

F  =  k  [x  -  x  ]    where  x  is  displacement. 
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The  velocity  is  the  time  derivative  of  the  displacement. 

From  the  harmonic  velocity  expressions  of  the  previous  section: 
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Equations  (l)  and  (2)  are  rewritten: 
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By  inspection  the  foir  pole  parameters  are 
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A-2.3   Damper  Four  Poles 
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The  performance  equations  are 

F  =  F 
1    2 

because  the  damper  possesses  no  inertia  and 

Fi  -  r<vi  -  V 

from  the  definition  of  damping. 
These  are  rewritten 


F  =  F 
1    2 


(1) 


(2) 


1   r  2    2 
By  inspection  the  four  poles  are 

eil=  ^  £12  =  °>  £21  =P  ^  622  =1' 
A-3  Parallel  Matrix  Combination 

If  n  systems  are  connected  in  parallel  then  a  single  four  pole 


(46) 

matrix  can  be  written  connecting  the  common  input  and  output  points.   If 
this  matrix  is  written  [S  ].  then 


-  B 


A  A£ 

B  B 

1  C_ 

B  B 


Let  the  r   matrix,  [S  ],  be  given  by: 


£11     £12 


r      r 
€21     €22 


where  r  varies  from  1  to  n.   Then,  in  (l)  above 


(1) 


(2) 


-£ 


~  I 


n 

r  m   1 


r  =  1 


(3) 


(4) 


(5) 


A  more  complete  discussion  is  given  in  reference  (6). 
A-k     Development  of  Four  Poles  of  General  System 
A-4.1   The  Machine 

As  stated  in  Section  2.4.1  the  machine  was  represented  by  a 
uniform  elastic  beam  driven  by  a  harmonic  force  generator.  This  system 
was  selected  because  it  was  mult i -re sonant  and  possessed  distributed 
parameters . 

The  four  poles  of  this  system  are  those  briefly  described  in 
Appendix  A-1.2  for  the  beam  in  the  preliminary  experiment.   These  four 
poles  were  e   =  i  e   =  — — ,  e        =  0,    and  e00  =  1.   From  the  method 


of  Reference  (k) ,   the  "beam  mobility  M(— )  will  be  derived. 
The  basic  equation  of  free  vibration  is: 

^x       ^  t 
where:  I  =  moment  of  inertia  of  beam's  cross  section 

E  =  modulus  of  elasticity 

jo  =  mass  density 

A  =  beam  cross-section  area 

y  =  transverse  displacement 

x  =  distance  along  length  of  beam 

L  =  length  of  beam 

t  =  time 
Assuming  a  harmonic  solution: 

y  -  f(x)eJu,t  (2) 

If  this  is  substituted,  equation  (l)  becomes: 

where 

and  \jJ   =  frequency  of  vibration. 
The  general  solution  of  (3)  is: 

f(x)  =  C,  (sinh^x  +  sin  ax)  +  C  (cosh^x  +  cos^x)  +  C,  (sinh 

ygx  -  sin^x)  +  Ci  (cosh^x  -t  cos^  x)  (h) 

The  boundary  conditions  for  a  free-free  beam  are: 

^2l-o  (5) 

dx 


d2f(L) 
dx2 


=  o  (6) 

4^=°  (7) 

dxJ 

4^=0  (8) 

Equations  (5)  and  (6)  require  no  shear  forces  at  the  ends  of  the  beam; 
equations  (7)  and  (8)  state  that  there  are  no  end  moments. 

When  the  boundary  conditions  are  applied  to  equation  (4),  the 
characteristic  frequency  equation  is  obtained. 

cos/6n  L  cosh/6n  L  =  1  (9) 

This  is  a  transcendental  equation.   Corresponding  to  each  eigenvalue  A 

will  be  a  value  of  f  (x).  The  term  f  (x)  represents  the  mode  shape  in 

which  the  beam  vibrates  at  the  frequency  defined  by  A    •  The  relationship 

between  f  (x)  and  A     is: 
n       rn 

fMx)  =  C2[cosh^en  x  +  cos^SQ  -o<n(sihh>Sn  x  +  sin^  x)  ]    (10) 

where 

sinh  A     L  +  sin  A     L 
.     _  rri r'n  ,   * 

jti  "  cosh  A     L  -  cos^   L  ^ 

^n        f^n. 

C  -   arbitrary  constant 


Define 


f  (x)  =  C  6    (x) 
nv      2  in 


Then  in  equation  (10): 

Q)    (x)  =  cosh  A     x  +  cos  A     x  ~eA    ( sinh  ^  x  +  sin  /3  x) 
in  r^n        rn      n      "n        r'n 

The  term  G)  (x)  will  satisfy  the  equation: 


(*9) 


0)(t)   dx  =  0    m^n 


f(|)^(x)^^4(L)]J 


dx 
It  is  a  normal  function  which  forms  an  orthogonal  set  in  the  interval 
O^xiL.  The  condition  of  orthogonality: 
L 

0 
However,  when  n  =  m: 

.L 

( 

0 

D.  Young  and  R.  Felgar  have  evaluated  these  (fy   functions  for  beams  with  a 
variety  of  boundary  conditions  (l6). 

With  the  values  of  (J)  and  the  corresponding  values  of  UJ  f    the 

forced  vibration  motion  is  considered. 

iu/t 
If  a  harmonic  force  F  e°   is  applied  to  the  beam,  the  displacement 

according  to  reference  (l)  is: 

FQeJwt     Y^X2n  (x) 
y(x)   =  — Z. ^5— — 2*+  rigid  body  motion  (12) 

T3      n  =  lu)  ~uJ 

n 

The  rigid  body  motion  is  f  (x)  eJ   .  The  velocity  is  the  derivative  of 
the  displacement: 

y(x)  =  — ^ Z_  -f ^+  Jwf0  (x)  e^  (13) 

"B       n  =  l(j„  -U>       U 
n 


(50) 
From  the  definition  of  mobility  at  the  driving  point  p: 

v(p)   y(p) 


=  iw  y       4n(l>)    .  ^0 


^  n  =  1  LU   -ui     0 

n 

As  a  rigid  body  the  beam  acts  like  a  mass  and  possesses  the  mobility  of 
this  lumped  parameter.   From  Appendix  A-2.1: 

-J 


roa-tois    l 

'lUUXJ-XUJ'      = 

h^u, 

Therefore,    in 

equation 

(lA)l 

JutfQ 

-J 

Fo 

mBa; 

The  final  expression  for  beam  mobility  is: 

m(p)=^  H  ¥^2-ib  (15) 

*ti  n  =  1  CJ2  -U2  ^ 

n 

Since  the  beam  in  the  machine  representation  was  considered  to  be 
driven  at  its  center,  the  four  pole  e   required  the  inverse  mobility 
at  the  center  of  the  beam.   This  was  expressed  from  equation  (15)  as: 

2    "B  n  =  lU2-u>2   V' 

n 

The  machine  could  now  be  represented  in  terms  of  its  four  poles. 
A-4.2   The  Mounting 

According  to  Section  2.4.2,  two  different  mountings  were  studied. 
The  first  was  the  standard  Portsmouth  mounting  which  was  analyzed  following 
the  method  of  Wright  and  Hagg  (17) • 

The  empirical  results  of  tests  on  actual  mountings  were  used  by 
Wright  and  Hagg.  The  mount  characteristics  were  approximated  by  use  of  the 


(51) 
electrical  impedance  analogy,  one  used  extensively  in  vibration  analysis 

(8).  The  mounting  was  represented  by  a  shunt  impedance  consisting  of  real 

and  imaginary  parts.   In  the  language  of  the  analogy,  impedance  is  defined 

as  a  force  to  velocity  ratio. 

l-f  (i) 

In  this  case  input  impedance  was  considered  for  the  approximation. 
From  the  analogy  mechanical  damping  corresponds  to  electrical  resistance 
or  the  real  part  of  the  impedance.  Mass  and  spring-like  properties 
correspond  to  electrical  reactance  or  the  imaginary  part  of  the  impedance. 
Thus  the  approximation  included  all  the  Important  lumped  parameters . 

These  parameters,  damping  constant  R  and  spring  stiffness  K,  were 
established  to  be: 

R==2^fQ^  (2) 

kx2  Wf2(l  +  f-) 
0 

k  = 3 2-  (3) 

gd  +  r^) 
t 

where : 

f  =  rated  natural  frequency  of  the  mount  at  its  rated  load,  cps, 

empirically  determined. 
Q  =  maximum  transmissibility  at  f  ,  dimensionless,  empirically 

determined . 
f  =  frequency  of  the  first  high  frequency  transmission  peak,  cps, 

empirically  determined. 

W 

—  -  rated  mass  load  on  the  mount,  lbm,  empirically  determined. 

S 

f  -  driving  or  exciting  frequency,  cps . 
The  shunt  impedance  Z  is  therefore  given  by: 


(52) 


Z  =  R  + 


K 


w 


The  impedance  analogy  indicates  the  following  equations  hold  for 
a  shunt  impedance: 


F.   n  F 
in    out 


(5) 


-    out        .    ., 
Z  =  — 1    open  circuit 

in 

V.   =  V  ,  ,  short  circuit 
in    out' 

Equations  (5)  to  (7)  are  translated  into  foiir  pole  parameter 

notation  as  in  Section  2 .2  r 

F  =  F 
1    2 

i z.h 

V-  =  V     station  2  free 


(6) 
(7) 


(8) 

(9) 

(10) 


Therefore,  the  four  pole  parameters  for  this  standard  Portsmouth 
mounting  are 

€11  m  **  €12  "  °'  £21  =  I  ^  €22  =  1" 
Combining  equations  (2),  (3),  and  (h)   the  following  expression  for  e 
is  obtained. 


u% 


"21   K 


1+  JQ0 


«0  +  1 


(53) 

(11) 


where  co  is  the  forcing  frequency,  rad/sec . 

The  second  mounting  considered  was  treated  similarly  except  that 
different  expressions  were  used  for  R  and  K.   These  were  derived  by  A.  0. 
Sykes  (l4)j 

R  =  ^2 (&  (12) 

ui(l  -  r)  L 

K  =  0(|)  (13) 

where  j 

S  *  mounting  cross-sectional  area 

L  =  mounting  length 

G  =  mounting  elastic  modulus 

r  =  mounting  damping  factor 

\/j  m  forcing  frequency 

The  elastic  modulus  and  damping  factor  were  given  as  empirical 
results  for  a  synthetic  high  damping  rubber  called  Thiokol  R.  D.  by  J.  C. 
Snowden  ( 18 ) . 

Approximate  formulae  were  determined  for  G  and  r  from  the  empirical 
data: 

G  =  72.6  UJ°'4  (lh) 

r  =  0.1U0'325  (19) 

Thus  R  and  K  could  be  represented  as  functions  of  frequency  and  the  new  four 
pole  parameter  e   evaluated. 

21   Gf(f-)  (A2  +  1) 


(5*0 


where 


A  = 


2r 


1  -  r 
A*4.3   The  Foundation 

As  stated  in  Section  2.4.3,  the  foundation  was  considered  a  beam 
resting  on  two  columnar  structures .  The  two  columns  presented  two  possible, 
identical  paths  for  vibration  transmission.  Since  the  paths  were  identical, 
only  one  was  studied  explicitly. 

The  four  poles  of  the  foundation  structure  were  a  series  combination 
of  the  individual  four  poles  of  a  beam  and  column. 


'11 


'12 


1 

0 


1 

m[l7 


121     22. 
The  combined  element  four  poles  are 


M, 


M. 


'11 


=  1  + 


WTJ'    €12   WJJ'    €21 


=  Mc,  and  a^   =  1. 


The  term  M(L)  is  the  transfer  mobility  of  the  beam  relating  the 
driving  point  with  the  junction  of  the  beam  and  column  structure.   Initially, 
this  beam  was  considered  as  end-supported  for  the  purposes  of  deriving  its 
mobility.   Correspondence  with  the  Admiralty  Engineering  Laboratory,  where 
an  experiment  driving  an  end-supported  beam  was  being  conducted,  indicated 
that  for  vibration  analysis  the  end-supported  beam  behaved  more  like  a 
free-free  beam.  For  this  reason  the  foundation  beam  was  treated  much  as 
the  beam  in  Appendix  A-4.1.   The  expression  for  beam  mobility  is: 

hl  *■ l   ,2    .2      dllJ 


"b  n  =  1  UJ      -U 


n 


V 


where  the  <fy ' s  are  the  normal  functions  at  the  driven  point  (— )  and  the  beam 
column  junction  (L) .  Again  these  functions  were  taken  from  Reference  (l6). 


(55) 
From  Reference  (l)  the  expression  for  the  lateral  deflection  of 
a  column  being  driven  by  a  harmonic  force  F  e    is : 


y(x)  =  f^  +   Z 

=  1  / (a 

n 


where 


p=   mass  density 

L  =  column  length 

A  =  column  cross;. section  area 

The  column  mobility  can  be  derived  using  the  same  methods  used  for  the 

beam  in  Appendix  A-^-.l.   The  velocity  is: 

.    ,       <=3°  n-1 

ju/t 


.,    ,         ,   ,     JuPb        frFOe            )  2(-l)2 

y(x)    =  jWfned        +  ^r / o 


v  n 
Now  the  definitions  of  mobility  and  the  rigid  mass  mobility  are  applied. 
From  this  the  mobility  of  the  column  is  obtained: 

eO        n-1 


M  =  ZJ_+  AJd 


L o 


2 
c   r    n  =  1  (u)  -u)  ) 


n 
A~4.^   The  Hull 

As  stated  in  Section  2..h.h,   the  work  of  D.  V.  Wright  was  consulted 
again  for  a  workable  analytic  representation  of  the  ship's  hull  (7)« 

Wright  considered  the  hull  a  series  impedance,  which  he  represented 
as  a  result  of  analyzing  empirical  data  obtained  from  typical  hull  impedance. 
measurements  (11,  20).  As  in  the  case  of  the  mounting  the  impedance 
consisted  of  a  real  and  imaginary  part,  the  resistive  and  reactive  properties. 
These  were  represented  as  complicated  quantities  dependent  on  nine  separate 
lumped  parameters . 


(56) 


R  = 


R     +  Rn  f  +  ^  +  -^ 
s  1  f 


(1) 


X  =  X     + 


p         p 


2   (123_  .  2nf } 

jtW  f        K              p    VjtW  f        K      ' 
m  s  s  *        p p 


+ 


3193         2ltf    '    ,    +  R2   (193         2*f  * 

1  +  RP  ( *w"T  "  k    } 
p         P 


Z  =  R  +  jX 


(2) 


(3) 


where  R  ,  Rn ,  R^,  R  ,  K  ,  K  ,  W  ,  W  are  the  nine  lumped  parameters  and 
s   1   2   p   s   p   s'   p 

f  is  the  driving  frequency. 

The  nine  parameters  were  considered  circuit  elements  in  the 
impedance  analogy,  making  up  the  series  impedance  Z. 


V. 
in  — 

-»■ 

Z 

V   , 
out  — 

— »- 

i 

i 
F 

Ol 

at 

F 

in 

' 

f 

i 

1 

The  impedance  analogy  gives : 

F. 

z   in 

z  "  v. 

in 

V.   =  V   . 
in    out 

F.   =  F  . ■  open  circuit 
in    out 


w 

(5) 
(6) 


Equations  (k)   through  (6)  are  translated  into  four  pole  parameter 
notation  by  the  method  of  Section  2 .2 : 


vi  =  v2 


F1    -   F     station  2  blocked 


(7) 
(8) 


(  ^  ) 


(57) 

- — f  = z  =  r-  <9) 

(M12)f 

The  four  pole  parameters  therefore  become : 

€n  =  1  (10) 

G12  =  Z  (11) 

e21  =  0  (12) 

e22  =1  (13) 

A-4.5   The  Water 

As  mentioned  in  Section  2  A. 5  the  water,  considered  only  as  a 
resistance  to  the  propagation  of  acoustic  waves,  has  a  specific  resistance 
R*  dependent  on  the  water  density  and  propagation  velocity.   This  can  be 
written,  (2) : 

R*=f=/>c  (1) 

where 

R*  =  Specific  resistance 

R  -   Resistance 

A  =  Effective  radiating  area 

fi  -   density 

c  =  velocity  of  sound  in  water 

Treating  the  water  as  a  dashpot-type  resistance  (velocity 
proportional),  the  four  pole  parameters  $?ere  obtained  in  Appendix  A-2-3* 

(2) 

(3) 

w 

(5) 


€11 

~ 

1 

€12 

= 

0 

£21 

= 

1 
R 

£22 

= 

1 

(58) 
R  was  evaluated  by  two  methods .   In  one  an  effective  area  A  was 

assumed  constant  with  frequency. 

R  =  ^cA  (6) 

The  hull  was  considered  also  as  a  vibrating  diaphragm.   The  radius 

r  of  the  effective  vibrating  area  is  given  by: 

r=£  (7) 

where  w  is  the  radiating  frequency. 
Therefore : 

R=-^  (8) 

Substituting  (6)  or  (8)  into  (h)   gave  the  desired  four  pole  para- 
meters . 

A-5  System  Properties  and  Systematic  Changes 
A-5.1   General  Description 

Four  variations  of  the  general  system  of  Figure  A-3  were  analyzed 
in  this  study.   These  variations  were  made  in  the  relative  sizes  of  the 
machine  and  foundation  structures  to  minimize  any  scale  effect  in  the 
trends  of  vibration  transmission.  Within  each  of  these  systems,  element 
properties  were  varied.   In  addition,  isolation  elements  were  inserted  at 
various  locations  in  the  systems.   From  these  variations  and  insertions, 
trends  to  provide  a  guide  for  preferred  design  practice  in  the  reduction 
of  vibration  transmission  were  established. 
A-5 .2   System  a 

With  reference  to  Figure  A-3  the  properties  of  System  1  were : 

1 .  Machine 

Weight  of  beam,  W^  =  50. 0  lbf 
Jd 

Mass  of  beam,  m  =  .1295  lbm 

£5 
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FIGURE       A-3 
GENERAL    SYSTEM 
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(60) 
Length  of  beam,  L  =  60.0  in 
Thickness  of  beam,  t  =  .75  in 

2 .  Mounting 

Rated  natural  frequency  of  the  mount  at  rated  load, 

Wo 

f0=-=l6.0cps 

Maximum  transmission  at  f  ,  Q_  =  7*0 

Frequency  of  the  first  high  frequency  transmission  peak, 

f .  =  ~  =  1000  cps 
t   2tt 

Rated  load  of  mount,  W  =  50. 0  lbf 

3.  Foundation 


Weight  of  beam,  W„  =  30.0  lbf 


Mass  of  beam,  hl  =  .078  lbm 

Length  of  beam,  L  =  36. 0  in 
Thickness  of  beam,  t  =  -75  in 
h .      Column 

Circular  Cross-section 
Length  of  column,  L  =  6.0  in 


Mass  of  column,  m  =  .039  Ihm 


Weight  of  column,  W  =15-0  lbf 


5 .   Hull 

The  nine  parameters  are : 
R  a  50  lbf -sec/in 


R 


^  .01  lbf -sec/in 


R  =  5000  lbf -sec/in 


(61) 
X  =  500  lbf -sec/in 

W  =10.0  lbf 
s 

K  =  10  lbf/in 

s         ' 

W  =  15  lbf 
P 

K  =  2  x  10  lbf/in 
P 

6 .  Water 

R  =  37^600  lbf -sec/in    (constant  with  frequency) 
A-5-3   System  2 

With  reference  to  Figure  A-3  the  properties  of  System  2  were: 
Identical  to  System  1  except: 

Foundation 

W^  =  100  lbf 

m„  =  .259  l"bm 

Jd 
L  =  60.0  in 
t  =  1.5  in 
A-5.4   System  3 

With  reference  to  Figure  A-3  the  properties  of  System  3  were 
Identical  to  System  1  except : 
Foundation 
W^  =  200  lbf 

ra  =  .518  lbm 

L  =  60.0  in 
t  =  3.00  in 
A-5.5   System  h 

With  reference  to  Figure  A-3  the  properties  of  System  k   were 


(62) 
Identical  to  System  1  except: 
Machine 
W^  =  500  lbf 

m_  =  1.295  lbm 

L  =  60.0  in 
t  =  3-75  in 
Foundation 
W  =  300  lbf 

m^  =  .78  lbm 

a 

L  =  60.0  in 
t  =  2.25  in 
A-5'6   Element  Insertions 
A-5.6.I  Thiokol  R.  D.  Mounting 

In  the  systematic  study  the  Thiokol  R.  D.  mounting  was  either 
substituted  for  the  Portsmouth  mounting  or  placed  in  parallel  with  it.  The 
properties  of  the  Thiokol  R.  D.  mounting  were: 

Mounting  cross  sectional  area,  S  =  9 «^3  in 

Mounting  length,  L  =  4.0  in 
A-5.6.2  Foundation  Damping 

Viscous  foundation  damping  which  was  constant  with  frequency  was 
first  considered.   The  damper  was  inserted  between  the  foundation  and  the 
column.   The  property  of  the  damper  was: 

Systeml-1: 

■n  •  4.       -o  onn    lbf -fee  C 

Damping   constant,    R  =  390  — ; 

System  2 : 

R  B  39o  iMz^SL 

JP  in 
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System  3' 

R  =  1560  ^=^ 

in 

System  k: 

lbf-sec 


R  =  1760 


in 


Damping  at  the  same  position  with  Thiokol  R.  D.  as  the  damping 

material  was  next  considered.   The  dimensions  of  the  Thiokol  R.  D.  were: 

2 
Cross  sectional  area,  S  =20.0  in 

Thickness,  L  =  1.0  in 

The  thickness  of  the  damping  material  was  such  that  its  compliance 
property  could  be  neglected. 
A-5-7   Systematic  Changes 

The  systematic  changes  in  the  structural  foundation  path  are  best 
followed  by  an  examination  of  the  individual  cases . 

Case  1.   System  1  used.   This  was  the  basic  case  for  System  1  with 

which  the  cases  with  property  changes  could  be  compared. 

Case  2.   System  1  used.  Mount  damping  was  increased  by  decreasing 

Q0  to  2.5. 

Case  3-   System  1  used.  Mount  stiffness  was  increased  by  setting 

W  =  201.062  and  U  =  12,566.384. 

Case  h.      System  1  used.  Mount  stiffness  was  decreased  by  setting 

U  =  50.266  andU  =  3,1^1.596. 

Case  5»  System  1  used.   The  Thiokol  R.  D.  mounting  replaced  the 

Portsmouth  mounting. 

Case  6.   System  1  used.  The  hull  damping  was  increased  by  setting 

R  =  100,  Rn  =  .02,  and  R0  =  2000. 

s      '   1      '      2 

Case  7*   System  1  used.   The  mount  damping  was  increased  by  setting 


(6k) 
Qn  =  2.5,  and  the  mount  stiffness  was  decreased  by  settintJ  =  50.266 
and  UJ  =  3} 1^1 -596.  This  case  was  thus  a  combination  of  cases  2  and 

Case  8.   System  1  used.   The  Thiokol  R.  D.  mounting  was  used,  and 
the  hull  damping  was  increased  by  setting  R  =  100,  R  =  .02,  and 

S  -L 

R  ■  1000.   This  case  was  thus  a  combination  of  cases  5  and  6. 
Case  9 •   System  1  used.   Mount  damping  was  increased,  mount  stiffness 
decreased,  and  hull  damping  increased.   This  was  accomplished  by 
setting  Q  =  2.5,  UJ  =  50.266,  U    =  3.IUI.596,  R  =  1000,  R  =  1.0, 

\J  \J  U  o  J- 

and  R  =  20,000. 

Case  10.   System  1  used.   Viscous  foundation  damping  was  inserted 

in  the  system  with  R  =  390* 

Case  11.   System  1  used.   The  Thiokol  R.  D.  and  the  Portsmouth 

mounting  were  placed  in  parallel. 

Case  12 .   System  1  used.   The  dissipative  resistance  of  the  water 

was  varied  with  frequency.   See  equation  (8)  in  Appendix  A-k.5. 

Case  13 »  System  1  used.  The  mount  stiffness  was  decreased  by  l, 

setting  u)   =  50.266  and  U,    =  3> 1^1 -596,  and  viscous  foundation 
0  t 

damping  was  inserted  in  the  system  with  R  =  390.  This  case  was 

thus  a  combination  of  cases  k   and  10. 

Case  Ik.      System  1  used.   Foundation  damping  with  Thiokol  R.  D. 

was  inserted  in  the  system. 

Case  l^.   System  1  used.   The  Thiokol  R.  D.  and  the  Portsmouth 

mounting  were  placed  in  parallel.   The  stiffness  of  the  Portsmouth 

mounting  was  decreased  by  setting  ijJ    =  50.266  and  U    =   3>  1^-1  •  596. 

This  case  was  thus  a  combination  of  cases  k   and  11. 

Case  l6.   System  1  used.  Hull  damping  was  increased  by  setting 

R  =  100,000,  R,  =  100,  and  R^  =  2  x  10  . 
s  '      1      '      2 
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Case  17.   System  1  used.  Foundation  damping  with  Thiokol  R.  D. 

was  inserted,  and  the  Thiokol  R.  D.  mounting  replaced  the 

Portsmouth  mounting.   This  case  was  thus  a  combination  of  cases 

5  and  Ik. 

Case  18.   System  1  used.   Foundation  damping  with  Thiokol  R.  D. 

was  inserted,  and  mount  stiffness  was  decreased  by  setting 

U    =  50.266   and  CJ  =  3, 1^1 -59^ .  This  case  was  thus  a  combination 

of  cases  k   and  14. 

Case  19 •   System  1  used.   Foundation  damping  with  Thiokol  R.  D. 

was  inserted,  mount  stiffness  was  decreased  by  setting  ul  =  50*266 

and  id     =  3>  1^-1  •  596,  and  mount  damping  was  increased  by  setting 

Q  =  2.5*   This  case  was  thus  a  combination  of  cases  19  and  2. 

Case  20.   System  2  used.  This  was  the  basic  case  for  System  2 

with  which  the  cases  with  property  changes  could  be  compared. 

Case  21.   System  2  used.  Mount  stiffness  was  decreased  by  setting 

WQ  =  50-206  and  bJ    =   3, 141.596. 

Case  22.  System  2  used.  The  Thiokol  R.  D.  mounting  replaced 

the  Portsmouth  mounting. 

Case  23 ♦   System  2  used.  Mount  damping  was  increased,  mount 

stiffness  decreased,  and  hull  damping  increased.  This  was 

accomplished  by  setting  Q  =  2.5,  Ui  =  50.266,  \jj=   3,141.596, 

R  =  1000,  Rn  =  1.0,  and  R„  =  20,000. 
s       *   1      '      2     ' 

Case  24.  System  2  used.  Viscous  foundation  damping  was  inserted 

in  the  system  with  R  =  390. 

Case  25.   System  2  used.   Foundation  damping  with  Thiokol  R.  D. 

was  inserted,  and  the  Thiokol  R.  D.  mounting  replaced  the  Portsmouth 

mounting. 

Case  26.   System  2  used.   Foundation  damping  with  Thiokol  R.  D. 
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was  inserted,  and  mount  stiffness  was  decreased  by  setting 

U£  =  50.266  and  u/  =  3,l4l.596. 

Case  30-   System  3  used.  This  was  the  basic  case  for  System  3 

with  which  the  cases  with  property  changes  could  be  compared. 

Case  31*   System  3  used.  Mount  stiffness  was  decreased  by  setting 

WQ  -  50.266  and  U  =  3,141.596. 

Case  32.   System  3  used.  The  Thiokol  R.  D.  mounting  replaced  the 

Portsmouth  mounting. 

Case  33-   System  3  used.  Mount  damping  was  increased,  mount 

stiffness  decreased,  and  hull  damping  increased.   This  was 

accomplished  by  setting  Q  =  2.5,  <J  =  50.266,  U    =  3,l4l.596, 

R  =  1000,  Rn  =  1.0,  and  R  =  20,000. 
s  1  2 

Case  3^»   System  3  used.  Viscous  foundation  damping  was  inserted 

in  the  system  with  R  =  I56O. 

Case  35 •   System  3  used.   Foundation  damping  with  Thiokol  R.  D. 

was  inserted,  and  the  Thiokol  R.  D.  mounting  replaced  the  Portsmouth 

mounting. 

Case  36 ♦  System  3  used.   Foundation  damping  with  Thiokol  R.  D. 

was  inserted,  and  mount  stiffness  was  decreased  by  setting 

L0Q  =  50.266  and  tjj    m  3,141.596. 

Case  50.   System  4  used.  This  was  the  basic  case  for  System  4 

with  which  the  cases  with  property  changes  could  be  compared. 

Case  51 »  System  4  used.  Mount  stiffness  was  decreased  by  setting 

U£  *  50.266  and  (J    =  3,141.596. 

Case  52.   System  4  used.  The  Thiokol  R.  D.  mounting  replaced  the 

Portsmouth  mounting. 

Case  53*   System  4  used.  Mount  damping  was  increased,  mount  stiffness 

decreased,  and  hull  damping  increased.   This  was  accomplished  by 
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setting  Q  =  2.5,  W=  50.266,  U  =  3,141.596,  R  =  1000,  R  =  1.0, 

and  R  =  20,000. 

Case  5^»   System  h   used.  Viscous  foundation  damping  was  inserted 

in  the  system  with  R  =  I76O. 

Case  55 »  System  k   used.  Foundation  damping  with  Thiokol  R.  fi. 

was  inserted,  and  the  Thiokol  R.  D.  mounting  replaced  the  Portsmouth 

mounting. 

Case  56.   System  k   used.   Foundation  damping  with  Thiokol  R.  D. 

was  inserted,  and  mount  stiffness  was  decreased  by  setting 

UL   =   50.266  and  U    =  3,1^1.596. 

Case  5$«   System  k   used.   Foundation  damping  with  Thiokol  R.  D. 

was  inserted,  mount  stiffness  was  decreased  by  setting  Ul  =  50.266 

and  W    =  3 j  1^-1  •  596,  and  mount  damping  was  increased  by  setting 

Q0  =  2-5. 
A-6  Sample  Calculations 

As  noted  in  Section  2.5«2,  the  computational  effort  involved  in 
the  thesis  was  so  great  as  to  render  it  impossible  without  the  use  of  a 
digital  computer.  However,  the  computational  effort  was  largely  needed 
in  iterative  calculations.  To  illustrate  the  complexity  and  length  of 
the  calculations  involved  in  solving  this  problem,  a  sample  calculation 
for  one  value  of  transmission  is  presented  here. 

For  example,  Case  1  properties  are  used  with  the  driving  frequency 
at  10.3  cps  or  65  rad/sec. 

The  unsimplified  four  pole  matrix  equation  for  Case  1  follows: 
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h 


M(fe 


Machine 


M 
Mounting 


M 


1  + 


M 


M(L)    M(L) 


Foundation 


H 


0 


Hull 


The  following  are  the  properties  for  this  case: 
1 .   Machine 


W„  =  50  lbf 


nL,  -  .1295  lbm 

L  -  60.0  in 
t  =  .75  in 
2 .   Mounting 

f  =  16  cps 


Q0  -  7-0 


ff  =  1000  cps 


W  =  50.0 
3.   Foundation 


W_,  =  30.0  lbf 
B   J 


m„  =  .078  lbm 
B 


L  =  36.O  in 
t  =  .75  in 


\     - 


Water 


(1) 


V, 
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h .  Column 

L  =  6.0  in 
m  =  .039  lbm 

W  =15.0  lbf 
c 

5-   Hull 

The  nine  parameters  are 

R  =  50  lbf -sec/in 

R  a  .01  lbf -sec/in 

R  a  5000  lbf -sec/in 

X  =  500  lbf -sec/in 

s  ' 

W  =10.0  lbf 
s 

K  =  10  lbf/ in 

W  =  15  lbf 
P 

K  =  2  x  10  lbf/in 
P  ' 

6 .  Water 

R  =  37^600  lbf -sec/in 
These  properties  are  also  listed  in  Appendix  A-5.2. 

The  individual  four  pole  parameters  must  be  evaluated  before 
combination  of  the  matrices. 
Machine : 

From  Appendix  A-U.l: 


M(f)    -  fF 


<t>n<l> 


"b  n  =  i  u2  -  ^    V 

n 


where 
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uj 


-*/ 


EI 


n         n       m 


<k  <t>  - 


B 


cosh, 


4° 

''  2 


L 


+  cos 


A: 


L 


sirih/2  L  +   sin^  L 
~n  ~n 


cosh>o  L  -  cos  A  L 
n  An 


sinh  - —  +  sin  - — 


A     and  Q)  (— )   are  listed  in  Reference   (l6).     uJ    is  calculated, 
n  Trr2  n 

^  =  3620 


u^  =  273 


W,  -  751 


W  =  5080 


(j^  =  -1.21565       ^6  =  ° 

d>2  =  0  cj>7  =  ~i.i+i4i2 


uj   =  1475        6JQ  =  6780  <|)    =  +1.42238        $Q  = 


W    =  2440         6U    = 


4 


9     r=     5 


V° 


(j)    =  +1.41412 


w5  =  3630        iuin  = 


10 


(|)    =  -1.41386       <t>     =  0 


10 


For  C<J=  65,  using  the  above  and  listed  properties 

10 i2,Ls 

'k\  -  JIM    ^  Tnjjj 

'O'  loot;       / o 


M( 


2y         .1295     ^1  ^2   _  ^  (65)(.1295) 

n 


M(f) 

=   - 

.1074 

J 

1 
MCg) 

=  9 

-30  j 

Mounting : 

From  Appendix  A-4.2: 


K  = 


(71) 
For  UJ=   65: 

(.1295X100.53)2  fi  +  fi^ 

K  =  i  ,  ?( 100 -537 
X  +  6283.2 

K  a  1285 

1_  _  (6^)(7)  rl  +  7Ji 
zM  "  1285    L^9  +  lJ 


|-  -  .00708  +  .o^96j 

ZM 


Foundation: 

From  Appendix  A-4.3: 

<='0    A  /L> 

M(L)  =  ^ 
m 


-> 

%^2J  %^ 

2          2 
1  UJ^    -UJ 

n 

B     n  = 

The  values  of  UJ  and  Q)  (L)  were  calculated  as  for  the  machine.   The  values 
n     Tn 


of  (j/(— )  are  the  same  as  for  the  machine. 


UJ±  =  758 

V 

10100 

4x 

=  +2.0 

§>6  =  92.0 

UJ2  =  2085 

6(J7  = 

li+080 

*2 

=  -2.0 

<t?  =  +2.0 

ui    =  4090 

H" 

18780 

fc 

=  +2.0 

fs  =  -2-0 

u)k  =  6770 

"9  = 

2^200 

<V 

=   -2.0 

^   =  +2.0 

u;   =  10100 

"V 

=   30100 

*5 

=  +2.0 

tl^O-fc.' 

For  uj  =  65  : 

10 

^WL^ 

•h  ^t 

.1 

M(L)  '  ^0W  2_  W2  .  ,225     (65K.078) 


M(L)   =  -   .2008  j 

W=M8 ' 
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oO         P 

M   =i^  >        2(-D  J 


c   m 
c 

'  n 


L =?  t   2    2N   m  u; 

n  =  1  (ur  -  u  )    c 
n 

From  Reference  (l): 


W  =  3-33  x  105 


U£  =  9.99  x  105 


W  =  16.65  x  105 


U^  =  23.31  x  105 

^c   =  29.97  x  105 
7 

For  6J  =  65  : 

5 

n-1 

Mc   ~  (.039)    f-T-i 

2(-D2 

(u£  -  ^225) 

3 
(•039)(65) 

Mc   =   -    .393   j 

Hull: 


From  Appendix  A-4.4 
Z  =  R  +  JX 


Rnaj       2itR_        R 
R  =  R     +     1     +_^+-P 


s       2ir  OJ  0  2 

p      p 

R„qn  +   (--01)  (65)    ,  2^(1000)        5000 

K   =    ^U    +  — +  £F^  +  2" 

R  =  146.704 
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tow        k         pW        k  ; 
x  =  x    +Z-1..--S.+ E P_ 


s         g 


1  +  R2(£ JiLl. 


P       *P 


6      (5000)2  (L,  -^ ^ 

x  =  500  +  ^ *T     I~~ 


(5000)  ((65)(l5) 


Water : 


2  x  10 

x  =  -  1^,898.3 

Z  «  146.7  -  14,898.3  j 

From  Appendix  A-4.5 
R  =  37,600 

|-  =  .0000266 
K 

Substituting  the  computed  four  pole  parameters  in  equation  (l) 
in  polar  form: 


~F1~ 

sa 

vl 

1      9.30 


0  1 


1      Ik 


Oh     1 


'Zf 


,902/f*~ 


.05/1. 


0 


428  1 


2.96/_0  4.98/f 

— 393/f       1 


2.66  x  10~5/o_         1 


(tt) 


This  equation  will  reduce  to  the  form: 


V. 


11     12 


£21    €22 


V, 


(3) 


With  reference  to  Figure  V: 


v6  =  o 


W 


From  the  four  poles  of  the  element  representing  the  water: 


(7<0 

F5  "  F6  (5) 

Substituting  equations  (k)   and  (5)  into  equation  (3): 

Fc    €11 
? 

In  section  2.5*1  the  transmission  T  was  defined  as: 
F, 


T  = 

20  log10  ^ 

Therefore : 

T  = 

-  -%o  ^ 

From  obtaining  equation  (3) J 
ei:L  =  5-51 

At  UJ=   65: 

T  *=  -  14,85 

A-7  The  IBM  709  Digital  Computer  Program 
A~7'l   Computer  Logic 

The  logical  procedure  comprising  the  program  written  to  solve  this 
problem  can  best  be  illustrated  by  the  computer  flow  diagram  in  Figure  A-4. 
A-7. 2   Scope 

The  cathode  ray  tube  oscilloscope  is  a  convenient  computer  feature 
which  produces  a  graphical  representation  of  the  problem  answers.   The 
computer  routine  used  to  program  this  feature  is  not  included  in  the 
program  listing.  Written) tin  FAP,  this  routine  projected  a  three  cycle 
semi-^og  grid  on  the  scope,  labeled  the  grid,  and  plotted  the  answers. 
Photographs,  35  mm*  in  size,  were  taken  of  the  scope  presentations.   These 
photographs  appear  as  results  in  Appendix  B. 
A-7 -3   Description  of  Subprograms 
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Figure  A-h 
Computer  Logical  Flow  Diagram 
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1.  MAIN  This  serves  as  the  control  program  calling  all  subprograms 
when  needed  in  the  solution. 

2 .  RDIN  This  subprogram  reads  in  the  data  including  the  system 
properties .  It  also  calls  the  routine  to  place  the  semi -log  grid 
on  the  scope. 

3-  MAT'l  This  program  computes  the  four  poles  for  the  first  element 
in  the  system.  In  this  case  the  four  pole  parameters  for  the 
machine  representation  were  obtained. 

k.     MAT  2  This  serves  the  same  purpose  as  MAT  1  for  the  second 
element  or  mounting  in  this  case . 

5 .  MAT  3  This  serves  the  same  purpose  as  MAT  1  for  the  third 
element  or  foundation  in  this  case. 

6.  MAT  h     This  serves  the  same  purpose  as  MAT  1  for  the  fourth 
element  or  hull  in  this  case. 

7.  MAT  5  This  serves  the  same  purpose  as  MAT  1  for  the  fifth 
element  or  water  in  this  case. 

8.  TRANS  This  program  takes  all  the  four  pole  parameters  computed 
by  the  various  MAT  routines  and  stores  them  for  use  in  the  matrix 
multiplication.   It  sets  up  the  matrices  prior  to  multiplication. 

9.  MAMUL  This  program  reduces  the  complex  matrices  to  one 
resultant  matrix. 

10.  ANSW  This  program  selects  the  matrix  elements  to  obtain  the 
answer  from  the  final  matrix  of  MAMUL.   It  also  provides  the  routine 
for  plotting  the  points  on  the  scope . 

11.  CHOMEG  This  program  tests  the  value  of  frequency  and  trans- 
mission to  insure  that  the  scope  points  are  numerous  enough  to  give 
an  adequate  plot.   It  also  changes  00  before  returning  to  the  start 
of  the  program. 
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12.  PRINO  This  program  prints  the  values  of  frequency  and 

transmission  to  conclude  the  program. 
A-7^   Program  Listing 

The  computer  program  listing  is  included  to  facilitate  its  possible 
future  use  in  design  investigations. 
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C      MAIN  RROukAM 

DlMtNMON  A(2b)»al2t>)*C<20)*D<20)»EUUU)»bU0U).ol  <  40U  )  » vM  2  0U)  »Vb 
1 (20)  »AA( 200)  »bb<  2 OU)  »CL I 20U)  »DU I 10U)  ♦ EL ( 10U) 

CALL  RUIN  ( A,b»L»u»AA»ub.CL»UU»tt ) 

G(  1  )  =  A( V) 

DO  210  L  =  lt<»00 
C      MATRIX  CALL  bEUUENCE 

CALL  MAU  (b»A.btb»U»L»AA»bb»LC»L>U»Et  ) 

CALL    MAT2     (  b  t  A  »b  tC  »D  »L  *  AA  ttiti  »CC  tDD  »tt  ) 

CALL    MAT3     (  b  »  A  #b  »C  »D  »L  »  AA  t  bb  »CC  »  0L>  *  tt ) 

CALL    MAl^f     (b.Atb»L*DtLtAA»bb»CC»UU»EE ) 

CALL  MA1S>  (b»A»b»L»U»L  . AA »bd »LL »UD » tt ) 
C       TRANSLATION  .StUUtNLt 

CALL  TRANb  (V»EE) 
C      MATRIX  MULTIPLICATION 

CALL  MAMUL  (  1  »J»K»V.A) 
C      ANSWERS 

CALL  ANbW  (Atbl»(j#t»LfV) 
C      CALL  FREQUENCY  CHANbt 

CALL  CHOMEG  (E»L»G»A) 

IF(G(L)-A<8> )  210»2<*0»2<f0 
210   CONTINUE. 
240   CALL  RRINU  (Gi»Gtfc»L) 

CALL  bXIT 

END 

TELE 


SUbROUTINE 

CALL 

FRAME 

CALL 

GRID 

CALL 

GRID 

CALL 

GRID 

CALL 

GRID 

CALL 

GRID 

CALL 

GRID 

CALL 

TITLE 

CALL 

VERT 

CALL 

HUR 

RETURN 

END 

SUbROUTINt  RDIN  (  A  »d  »C  »L>  »  AA  »bd  »LL  t  Du  »  tt  ) 

DIMENSION  A  <  2:>  >  »b  (<;:>)  »L  (20)  »U  I  20)  »AA(  <:UU  >  »bb<  <iOU)  *CC  I  200)  »UD<  lUU  )  * 
1EE(100) 

READ  213»((A(I)»l=l*lU)*(b(I)*I=l*2O)»(C(I)tI=l»20)»(D(I)»I=l»^O)t 
1(DD( 1  )  t  1  =  1.10)  ) 
213   FORMAT  OE13./) 

CALL  TELt 

RETURN 

END 

SUbROUUNt  MAfl  Ib.A.d.C.U.L.AA.ob.LLtDU.tt ) 
btAM  MATRIX 

OIHtNilON  bl^UO )  iMl<:3)  *U«  a'j)   »  L  (  <_  U  )  »  L>  I  <:  U  )  *AA«^UO)  »bdl 2UU )  »Cll<iUO)  * 
1DD( 100) »tE( 100) »Vb(20) 
VV3=1.0/ (b(2  )»b(l.  ) » 


(7*0 

DO  20b  I  =  i»  10 

IF(D( 1 )**2-0(L)**2 >  201»202»20l 

201  VB<  1  )=G(L)*C<  I  1**2/  (U(  (  )**2~0<L)*  'f2  ) 
GO  TO  20b 

202  IF  (G(U*C( 1 )**2-0«0)  204»203»204 

203  VB(I)=0.0 
GO  TO  203 

204  VB( I )=A<  /  ) 
20t>       CONTlNUt 

V78*Vb(i )+vb«2»+Vb(3)+Vb(4)+Vb(b)+Vb(6)+vtt( n+Vb(a)+vb( v)+vb( iu) 
v/y=v7«/bu ) 

V229=V79~V43 
IF(V229-0.0)  <iOb»20d»20  / 

206  V329=A{b) 
GO  TO  209 

207  V329=A(2) 
GO  TO  209 

208  V3^V=0»0 
V32d=A(  1  ) 
GO  TO  212 

209  V32B=1.0/AbSF< V22V ) 
212   EE(1)  =  \Z52H 

EE(2)  =  V329 

RETURN 

END 

:      THIS  IS  MOUNT 

SUBROUTINE  MAT2  ( G» A »b »C »D .L »AA »bB »CC »DD »EE ) 

DIMENSION  GUOO)  ♦A(2b)tb(2b)»C<20)»D<20)»AA(200)  »bB(200)  ♦CC(200) » 
1DD( 100) »EE(100) 
VA  =  <b(4)**2*B(b)  )*(  (Bib  )+b(L)*A(l ) ) / ( b I  6 ) +b (4 ) *A < 1) )  ) 
RE  =  (G(L)*B(3)  )/(VA*(B<3)**2  +  1.0) ) 
VE  =  RE*B(3) 

EE(2)  =  SQRTF( VE**2+RE**2) 
EE(4)  =  ATANF(VE/RE) 
RETURN 
END 

SUBROU TINE    MAT2(b»A»otC*D»L»AA»dd.CC>DD»EE) 

DIMENSION    GUOO)  *a(  cb)  »b(2b)  »C(2U)  »U(20)  »AA(200)  »BB(200)  »CC(200) » 
1DD( 100) »EE(iOO) 
GEE    =    72.t>*EXPh  CD.4*L0GF  (  G  {  L)  )  ) 
RAR    -    0.1*EXPF(0.32b*LOuMb(L  )  )  ) 
V43    =    RAR**2 
IF(V43-1.0)     410»411»410 

410  RE    =     (B< 19)/B<20)  )  <HbEE/G(L) )*(2.0*RAR/(1.0-V43) ) 
VA    =     (B( 19)/B(20) )*(GtE/G(L) ) 

EE(3)     =     1.0/SORTF(Rt*#2+VA**2 ) 
EEU)     =    ATANF(VA/RE) 
GO    TO    412 

411  EE(3)     =    0.0 
EE<4)     =    0.0 

412  RETURN 
END 


410 


411 
412 


(80) 

SUbROUl  INt  MA  I  /  (G.AtB.C  »U » I »AA »bb »CC »DD » EE ) 

DIMENSION  GUOO)  .A(23) .B(23) »C(20) »D(20) »AA(200)  .BB(200) »CC(200)  . 
1D0( 100) »EE ( 10U) 
OEE  =»  /2»0'it:APl-  (U.4*L0bh  ( <j  (  L)  )  ) 
KAK  =  0»1*EXPI-  (0»32b*LOljF  ( b  (  L)  )  J 
V4  3  =  KAU**2 
lF'(V'+3-l«0)  4  1U  »'>!  1  itiU 

RU  =  (b(  19  )/B<  <iO)  )  *l(jtt/o(L>  )*(2.0*RAR/<  1.0-V43)  ) 
AA  =  IBIltfl/HI'iUDMbtL/OlLlI 
GO  TO  412 
RR  =  A  (7) 

VA  =  (BU)**2*l3(bM*((B(6)+G(U*A(l))/(B(6)*BU)*A(im 
RE  =  (G(L)*B(3  )  )/  <VAMB(3)**2  +  1.0)  ) 
VE  ■  RE*B(3) 

V47  =  SUk  TK  IkK'^+AA*^  ) 
V48  -  ATANF (AA/RRJ 
V49  =  SQRTF (RE#*2+VE**2 ) 
V30  =  ATANF (VE/RE ) 
TOPI  =  V4/»V49 
T0P2  =  V4b  +  Vi>0 

B0T1  =  SuRIF  {  (l<EH<k)**2+(  AA+Vfc)  **2) 
B0T2  ■  Al AN*  ( ( AA+VE )/<RE+RR) ) 
EE(3)  =  T0P1/B01 I 
EE(4)  =  T0P2-B0I,' 
RETURN 
END 


C      THIS  IS  PUUNUAI 1UN 
C       MATRIX  THKtL 

SUBROUTINE  MA  I  J  ( (j»  A  tB  iC  »D  tL  t  AA»BB  t  Ct  •  DD.fc.E  ) 

DIMENSION  Gl4OU)«A(2i>)»Bl2i>)»C(2U)»U(2U>tAA(2UO)»BB(20O)tCC(2U0)» 
100(100) »EE( 100) »VB(^0) 

V4  3  =  l»0/(b(7)  *<j<L)  ) 

DO  303  I  =  1»10 

IF(D{ I+i0)**2-G<L)**2)  301  »  302.301 

301  VBII+10)  =  (G(L)«U 1 ) «C( 1+10) )/(0{ I+10)**2-G(U**2) 
GO  TO  303 

302  IF  (G(L)*C(I)*C(  I+.JLO)-0«0)  30<+»303»304 

303  VB(I+10)  =  0.0 
GO  TO  303 

304  VB< 1+10)  =  A( 71 
303   CONTINUE 

V76  =  Vb(ll)+Vb(i2)+VM( 13)+Vb( 14)+Vb( 1 3 ) +Vb ( 16 ) +Vb ( 17)+Vb( 18)  + 
1VB( 19 )+VB(20) 
V79  =  V78/b(7) 
V229  =  V79-V43 
IF(V229-0.0)  306,308. 40  7 

306  EE(6)  =  A(6) 
GO  TO  309 

307  EE(6)  =  A(^) 
GO  TO  309 

308  EE(6)  =  0.0 
EEI5)  =  A( I) 
GO  TO  312 

309  EE(3)  =  l.U/Abol- (V2^"y  ) 


(81) 


312 


RETURN 
END 


301 

302 
303 

304 
305 


306 
307 
308 


309 
312 


301 

302 
303 

30<+ 

30b 


THIS  lb  HUH  l-UUNUAl  ION  WIIH  DAMRtR  (-OLtOWlNij 

MATRIX  FHKhE 

SUBROU1  INfc  MA  I  J  (G»A»M»C»D  »L  »AA  »BB  |U  »DL>  •  tt  ) 


DIMENSION  (j(  ^00  )  »A< 
1DDU00) tEEUOO) «VB(2U) 
V43  =  1.0/ (B(  7) *GU  »  ) 
DO  30b  1  =  1»10 
1F(D(  1  +  10)  **2-U(U  "  ^) 


)  .hi  /',  )  .(J^ii)  »D(  ^U  )  »AA(  <^00)  »!JB(  2  00  )  »CC  (200)  » 


*01»302»301 
(G(t)*U  I  )  "<•  (  1  +  1  0)  >/(u(  1  +  10)  **2-G(L  >**<>) 


O.o)  304»303t304 


VBU  +  10)  = 
GO  TO  30b 

IF  (ti(L)*l( 1 ) »K 1+10) 
VBU  +  10)  =■  0.0 
GO  TO  30b 
VBU  +  10)  =  A(  7) 
CONTINUE 

V78  *  VB(  11  )+VB(  12  )+VH<  1  3  )+\/H(  1<0+VB<  lb)+Vb(  ib)+VB(  1  /  )+VB(  18)  + 
1VBU9)+VB(20) 
V79  =  V7b/B<  U 
V229  =  V/9-V43 
IF1V229-0.0)  -JObtHObt-^U/ 
EE(6)  =  A(6) 

309 

=  A(2) 

309 

=  0.0 

=  A  (  7  ) 

312 

=  1.0/AtHol-  (  V^y) 
l.O/Al 10) 
l.U+lhthl-ttliol) 


GO  TO 
EE(6) 
GO  TO 
EE(6) 
EE(b) 
GO  TO 
I  i  (  5  ) 

(16)  = 
tt(14)  = 


VA  =  (EE(b  )  ;;tt  I  lb  )  )  -Mm  lttlt>() 

RE  =  1.0+(tt (b  )*tt ( lb)  )  •  'LOSh  ( it (to  )  ) 

EE(  lb)  =  ATANI-  (VA/Rt  ) 

RETURN 

END 

FOUNDATION  WITH  UAMFlR  VmKIAULl  Wl  In  IktUUtNCY 
MATRIX  THRtt 

SUBROOIINt  MA  l  3  ('o»A»Bt'  »0  »L  •  AA  ♦Bt)  t  V.L  »DD  » tt ) 

DIMENSION  G1400 ) »A(2b ) »B(^b) »t (20) »D( 20) »AA(200) »BB(200) »CC(200) « 
1DD(  100)  »Et (  100)  »VB(^0) 
V4  3  =  1.0/(B<7)«o(D) 
DO  30b  I  =  1,10 

IF(DU  +  10)**2-ML>#*2)  301,302t301 

VBU  +  10)  =  (G!L)*CU  )*CU  +  10)  )/(DU+10)**2-<9<L)**2) 
GO  TO  30b 

IF  (G(L)*CU  )*CU  +  10)-O*0)  304»303i304 
VB(  I  +-10  )  =  0.0 
GO  TO  30b 
VB(  1  +  10)  =  A(  /  ) 
CONTINUE 
V78  =  VB(  11  )+VB(  i2  )+Vd(  13  )+Vh(  m)+VB(  )  b  )  +  VB  (  lo  )  +  VB  (  1  f)  +  VB  (  lti  )  + 


(82) 


500 


501 
C02 


lVb(  19)+Vii<20) 
V79    3    V  /b/b(  /  ) 
V229    =    V79-V4  3 
IFJV229-0.0)     306*30«»307 

a  Alb  ) 

30V 

=  A<2  ) 

309 

=  0.0 

»  A(  7) 

312 

=  1.0/ABbl-  (  V229) 

72.6*tXf^  <0.4»Luoe  (b(L) ) ) 

0»1*EXPF (0.32  5"LUbh IblLIII 

KAR**2 
IF (V44-1.0)  500*501*500 

RE  =  b(l  )  /A(  10)*(GEE/G(L)  )*(2.0*RAR/(1.0-V44)  ) 
EE(16)  =  1.0/Kt 
GO  TO  502 
EE(16)  =  0.0 

EEI14)  =  1.0+<EE(5)*EE( 16) ) 
VA  =  (EE(5)*EE( 16)  )*SIM  (EE<6>  ) 
RR  =  1.0+(EE(5)**EE(16) )*COSF(EE(6) ) 
EE(15)  =  ATANFfVA/RR) 
RETURN 
END 


306 

EE(6) 

GO  TO 

307 

EE(6) 

GO  TO 

308 

EE<6) 

EE(5) 

GO  TO 

309 

EE<5) 

312 

GEt  = 

RAR  = 

V44  = 

401 

402 
405 


406 
407 


THIS  IS  COLUMN  HULL  MATRIX  4 

SUBROUTINE  MAT4  ( G ♦ A . b *L »D • L »AA »db » CC »UO tEE ) 

DIMEN6I0N  G(4U0)tA(25)»b(i'5)»C(20).D(20)»AA(200)  »bB(200)  »CC<200>  ♦ 
1D0(  100)  *EE(  100)  *Vb<<d0) 

V4  3  =  liO/(d(8)*(i(UI 

DO  405  I  =  1»5 

IF  (DD(I )**2-G(L)**2)  401*402*401 

VB(I)  =  (G(L)*2.0*DD( 1+5) > / ( B ( 8 )* ( DO ( I ) **2-G ( L > »*2 ) ) 

GO  TO  405 

VB(  I  )  =  A(7) 

CONTINUE 

V78  =  Vb( 1  )+VB(3)+VB(5)+V43 

V79  ■*  VB(2)+VB(4) 

EE(7)  =  SQRTF(V78**2+V79**2) 

IF (V79-0.0)  407*406 .407 

EE(8)  =  A(6) 

EE(8)  =  ATANF(V7b/V79) 

V45  =  B(9)  +(B( 11 )*G(L) ) / A ( 3 ) +A ( 3 ) *B ( 12)/G(L) 

V46  =  386.0/B( 17)*G(L)-o(L )/B( 16) 

V45  =  V45+b( 15)/< 1. 0  +  b( 15  )'  2  V46-  >2) 

V4  7  =  b( 10)+b< 14) *GIL)/3B6.0-B( 13)/G(L)+(B( 15 ) **2*V46 ) / ( 1.Q+BU5) 
1*2*V46**2) 

EE(9)  =  SORTF( V45 **2+V47 * « 2 ) 

EE(10)  =»  AlANh"  (V47/V45  ) 

V50  =  EE(7)*EE(9) 

V51  =  EE(8)+EE(  10) 

V52  =  V50*CUSF (Vbl) 

V53  =  V50*SINF(V51) 


(83) 

V52  =  1.0+V52 

EE(ll)  =  SQRTF(V52**2+V53**2) 

EE(12)  *  ATANF<V53/V52) 

RETURN 

END 

C      MATRIX  FIVE 
C      THIS  IS  WATER 

SUBROUTINE  MAT5  (  Gt  A  »Bt  C  »D  »L  »AA  »bb  »CC  »0UiEE  ) 

DIMENSION  0(400)  »A(25)»B(25)tC(20)»Ol2O).AA(20O)  »BB  (  200  )  #CC  (  200  )  » 
1DD( lOO).EE(lOO) 

EE(13)  =  1.0/BU8) 

RETURN 

END 

■   SUBROUTINE  MAT5  < G »A »B »C »D »L »AA »bB »CC »DD »EE  ) 

DIMENSION  GUOO)  *A(25)»B(25)»C(20)»D<20)»AA(200)  »BB(200)  »CC<200)  . 
lDD(lOO) .EE(IOO) 
EEU3)  =  G(L)**2/B(1) 
RETURN 
END 

C      MATRIX  MULTIPLICATIUN 

SUBROUTINE  MAMUL  (ItJtktVtA) 
DIMENSION  V<200)»A(25) 
1  =  1 
J  =  9 
K  =  68 
101   K=K+2 

V41=V( I >*V( J) 

V42=V( I+l)+V( J+l) 

V43=SINF(V42)*V41 

V44=C0SF(V42J*V41 

V45»V( I+2)*V< J+4) 

V46=V< I+3)+V( J  +  5) 

V47=SINF( V46) *V45 

V48=C0SF(V46)*V45 

V49=V44+V48 

V50=V43+V47 

V(KJ=SQRTF(V50**2+V49**2) 

IF(V44+V48)  50»25»50 

50  V(K+1)  *  ATANF(V50/V49) 
GO  TO  26 

25  IF1V43+V47)  51»  27t  51 

51  V(K+1)  *  A(6) 
GO  TO  26 

27    V(K+1)  =  0.0 

26  IF(K-72)  1099»  102.  1100 

1100  IF(K-74)  1099»  3»  1101 

1101  IF(K-76)  1099»4»1102 

1102  IFU-78)  1099»5»1103 

1103  IF(K-80)  1099#6»1104 

1104  IFU-82)  1099»7.1105 

1105  IF(K-84)  1099»8»1106 

1106  IF(K-86>  1099t9»1107 


(84) 

1107  IF(K-88)  1099,10,1108 

1108  IF(K-90)  1099,11.1109 

1109  IF(K-92)  1099,12,1110 

1110  lF(K-9<f)  1099»13»1111 

1111  IF(K-96)  1099,14,1112 

1112  IF(K-98)  1099,15,1113 

1113  IF(K-IOO)  1099,16.1099 
1099  I  =  1+16 

J  *    J+16 
GO  TO  101 
102   V102  =  (V(70)*V(72) ) *SINF < V ( 71 ) +V{ 73  ) ) 
V103  »  (V(70)*V(72) )*  COSF(V(71)+V(73)  ) 
I  =  1-16 
J  =  J-14 
60  TO  101 

3  I  =  1+20 
J  =  J+14 
GO  TO  101 

4  V104  =  (V(74)*V(76) )*  S INF ( V ( 75 ) +V ( 77  )) 
V105  =  (V(74)*V(76> )*C0SF<V(7i  )+V(77) ) 
V(106)  =  V104+V102 

V(107)  =  V105+V103 
I  *  1-16 
J  =  J-16 
GO  TO  101 

5  I  =  1+12 
J  =  J+16 
GO  TO  101 

6  V102  =  (V(78)*V(80) ) *S I NF < V ( 79 > +V ( 8 1  ) ) 
V103  =  <V(78)»V(80) )*C0SF(V(79)+V(81> ) 
I  =  1-12 

J  =  J-14 
GO  TO  101 

7  I  =  1+16 
J  =  J+14 
GO  TO  101 

8  V104  =  (V(82)*V<84)  )  *S  INF ( V ( 83 ) +V ( 85 ) ) 
V105  =  (V(82>*V(8M ) *COSF ( V ( 83 ) +V < 85 ) ) 
V<  108)  =  V104  +V102 

VI109)  =  V105+V103 
I  =  1-20 
J  =  J-16 
GO  TO  101 

9  I  =  1+16 
J  =  J+18 
GO  TO  101 

10  V102  =  (V(«6)*V(88) ) *b 1 NF (VU7)+V(U9 ) ) 
V103  =  IV(86)*V(88) )*COSF ( V< 87 )+V ( 89 ) ) 
I  =  1-16 

J  =  J-16 
GO  TO  101 

11  I  =  1+20 
J  =  J+16 
GO  TO  101 

12  V104  =  (V(y0)*V<92))*i>INMV(91)+V(93)) 


(85) 

V105  =  (V(90)#V(92)  )*COSF(V(91)+V(93)  ) 

V(110)  =  VI 04  +  y/102 

V(lil)  =  V10'^+  V103 

I  =  1-16 

J  =  J-ld 

GO  TO  101 

13  I  =  1+12 
J  ■  J  +  18 
GO  TO  101 

14  VI 02  =  (V(94)*V<96) ) *S1NF ( V( 95 ) +V < 97 ) ) 
V103  =  (V(94)*V(96) ) *COSF ( V ( 95 )+V(97 ) I 
I  *  I-  12 

J  ■  J-16 
GO  TO  101 

15  1  =  1+16 
J  ■  J+16 
GO  TO  101 

16  V104  ■  (V(98)*V(100) )*SINF ( V< 99 )+V( 101 ) ) 
V105  =  (V(98)  *V(  100)  >  «COSF(V(99)+V<101)  ) 
V(112)  =  V104  +  V102 

V(113)  =  Vl0b+V103 
I  =  1+85 
J=   J  +  6 
K  =  ^+20 

17  K,  =  K.  +  2 

V43  =  SQRTF(V<  I  l**i!+Vli+l)**2) 
IF(V(  1  +  1  )+0.0)  5>2»2d»t><i 
i?2    V44  =  ATANMV(1»/V(1+1I) 
GO  TO  29 

28  lF(V(I)+0.0)  53»30.b3 
53    V44  =  A(6) 

GO  TO  29 
30    V44  =  0.0 

29  V(K)  o  V43*V( J)*S1NF(V44+V< J+l) ) 
V(K  +  1)  a  V43*V( J)*C0SF(V44+V( J+l)  ) 
1F(K-124)1098. 18*1115 

1115  IF(K-126)109B»19»1116 

1116  IFU-128)  1098»20»111  7 

1117  IFU-130)  1098»21#1118 

1118  IF(K-132) 10Vd,22»liiy 

1119  IFU-134)  10V8.23»1120 

1120  IF(K-136) 1098*24. 10^8 
1098  I  =  1+4 

J  =  J+4 
GO  TO  17 

18  V138  =  V(122)+V(124) 
V139  =  V( 123)+V(125) 
I  =  1-2 

J  =  J-4 
GO  TO  17 

19  1  =  1+4 
J  ■  J+4 
GO  TO  17 

20  V140  =  V( 126)+V( 128) 
V141  =  V( 127)+V(  129) 


(86) 


I  =  I-S 
J  =  J-2 
GO    TO    17 

21  1=1+4 
J  =  J+4 

GO    TO    17 

22  V142=V(130)+V(132) 
Vl43=V(131)+V(133) 
1  =  1-2 

J  =  J-4 
GO    TO    17 

23  1=1+4 
J  =  J+4 

GO    TO    17 

24  V144  =  VU34)+V(136> 
V145=V(135 >+V( 137) 
RETURN 

END 

THIS  IS  TRANSLATION 

SUBROUTINE  TRANS  (V»EE) 

DIMENSION  V(200) »EE( 100) 

V(l)  =  1*0 

V(2)  =  0.0 

V(3)  =  EEU) 

V(4)  =  EE(2) 

V<5>  ■  0.0 

V<6)  =0.0 

V(7)  =  1.0 

V(8)  =  0*0 

V<9)  =  1.0 

V(10)  =  0.0 

V(ll)  «  0.0 

V(12)  =  0.0 

V(13)  =  EE(3) 

V(14)  =  EE(4) 

V(15)  =  1.0 

V( 16)  =  0.0 

V(17)  =  1.0 

V(16)  «  0.0 

VU9)  =  EE(5) 

V(20)  =  EE(6) 

V(21)  =  0.0 

V(22)    =  0.0 

V(23)  =  1.0 

V(24)  =  0.0 

V(25)  =  1.0 

V(26)  =  0.0 

V(27)  =  EE(9) 

V<28)  =  EE(10) 

V(29)  =  EE(7) 

V(30)  =  EE(8) 

V(31)  «  EE(ll)  ' 

V(32)  =  EE(12) 

V(33)  «  1.0 


(87) 


V(3<t)  « 

0*0 

V(35)  » 

0.0 

V(36)  » 

0*0 

VI37I  = 

EE(  13  » 

V(38)  = 

0.0 

V(39)  * 

1.0 

VUO)  ■ 

0.0 

RETURN 

END 

*      FORTRAN 

C      THIS  IS  TRANSLATION 

C      THIS  IS  FOR  FOUNDATION  WITH  DAMPER  FOLLOWING 

SUBROUTINE  TRANS  (V»tt) 

DIMENSION  VI200) #EE( 100) 

V<  1  )  «  1.0 

V(2)  »  0.0 

V<3)  »  EE (  1  > 

VU)  =  EE(2) 

V(5)  a  0.0 

V<6)  «0.0 

V(7)  =  1.0 

V(8)  «  0.0 

V(9)  »  1.0 

VUO)  =  0.0 

V(ll)  =  0.0 

V(12)  =  0.0 

V(13>  =  EEI3) 

V(14)  =  EE(4> 

V(15)  =  1.0 

V(  16)  =  0.0 

V(17)  =  EE(14) 

V(18)  =  EE(15) 

V(19)  =  EEI5) 

V(20)  *  EE(6) 

V(21)  =  EE(16) 

V(22)  =  0.0 

V(23)  =  1.0 

V(24)  =  0.0 

V(25)  =  1.0 

V(26)  =  0.0 

V(27)  =  EE(9) 

V(28>  =  EE(IO) 

V(29>  *  EE(7) 

V(30)  =  EE(8) 

V(31)  =  EE(ll) 

V(32)  =  EE(12) 

V(33)  =  1.0 

V(34)  =  0.0 

V(35)  =»  0.0 

V(36)  =  0.0 

V(37)  =  EE( 13) 

V(38)  =  0.0 

V(39)  *  1.0 

V(40)  =  0.0 
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RETURN 
END 

SUBROUTINE  ANSW  ( A»G1 »G»E »L tV) 

DIMENSION  A (25) »G1(400) tGUOO) »E(400) »V(200) 

V138  =  V( 122)+V( 124) 
V139  =»  V(  123)+V(125) 

E(L)  =  (A(5)/AI4) )*LOGF(1.0/SQRTF(V138**2+V139**2) ) 
GKL)  =  G(L)/A(3) 
G2  =  LOGF(GKL)  )/A(4) 

MFREO  =  ( ( (G2-1. 01*920. 0)/3.0)+102«0 
MT  =  (( (E(L)+70.0)*870.0)/140.0)+152.0 
CALL  BRlPIC(MFREOtMT) 
CALL  BRIP1C  (MFREO»MT) 
RETURN 
END 

SUBROUTINE  CHOMEG ( E »L tG »A ) 
DIMENSION  EUOO)  »G<400)  »A(25) 
IF(G(L)-A(9) )251»250»251 

251  V502  =  ABSF(E(L)-E(L-1) ) 

V1012  =  ABSF(LOGF(G(L) ) -LOGF ( G( L-l ) ) ) 
IF(V502-5«0)254»253»25  3 

254  IFIV502-2.5  )250»255»255 

255  IF(V1012-0.02878)250»253»250 
253   G(L+1)  =  EXPF(LOGF(G(L) )+. 02878) 

GO  TO  252 
250   G(L  +  1)  =  EXPF(LOGF(G(L)  )  +  . 05757) 

252  RETURN 
END 

SUBROUTINE  PRINO  (G1»G»E»L) 
DIMENSION  GK400)  »G(400)  »E(400) 
PRINT  220.  (G( I ) »G1( I ) »E( I ) »I*ltL> 
220   FORMAT(8H  OMEG  =  F8.1»9H   FREO  =  F7.1.6H   T  =  E13.6) 
CALL  FRAME 
RETURN 
END 
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APPENDIX  B 

PRESENTATION  OF  RESULTS 

The  following  photo  plots  were  produced  by  the  Cathode  Ray  Tube 
Recorder  of  the  IBM  709  Digital  Computer.   The  computed  points  were 
connected  by  solid  lines  by  the  authors  to  clarify  the  plots . 

Due  to  the  space  limitations  of  the  scope  presentation  the 
abscissa  or  frequency  scale  is  scaled  in  integers  for  three  decades. 
The  first  decade  on  all  of  the  plots  extends  from  10  to  100  cps,  the 
second  decade  from  100  to  1000  cps,  and  the  third  from  1000  to  10,000 
cps .   The  ordinate  presents  transmission  in  decibels  varying  from 
-70  to  +60. 
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